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1  Executive  Summary 

The  aim  of  this  section  is  to  provide  an  overall  sumnoaiy  of  the  research  carried  out  with  the  full 
or  partial  support  of  the  contract  Although  formal  support  was  only  provided  for  one  PhD 
research  student  A  J  Sant  in  practice  other  students,  notably  M-J  Kim,  N  C  Bruce  and  D  N 
Qu  all  beneritted  indirectly  from  ERO  support  and  contribute  to  the  project 

Background 

To  set  the  research  in  context  it  is  worthwhile  reviewing  the  work  of  ourselves  and  others 
prior  to  mid- 1987,  concentrating  on  enhanced  backscattering  firtxn  randomly  rough  surfaces. 
There  is,  of  course,  a  vast  literature  on  the  scattering  of  electromagnetic  waves  by  randcnnly 
rough  surfaces.  Prior  to  1987,  the  vast  majmity  of  theoretical  studies  were  performed  using 
analytical  methods  and  it  was  therefore  necessary  to  make  approximations  in  order  to  get  any 
kind  of  answer  at  all  for  the  scattering  cross-section  (ensemble  average  intensity).  In  the 
great  majority  of  treatments  assumed  either  a  perturbation  situation  (for  which  die  rms  surface 
height  a  typically  has  to  be  less  than  X/100,  where  X  is  the  wavelength)  or  one  where  single 
scattering  and  the  Kirchhoff  approximation  applied  (i.e.  surfaces  for  which  a/x«l,  where  t  is 
the  correlation  length  of  the  surface  roughness)  —  this  is  also  called  "Beckmann  theory"  or  the 
"physical  optics"  approach.  The  first  pi^ction  of  enhanced  backscattering  was  made  by 
McGum  et  al^  using  an  analytical  approach  in  the  perturbation  situation;  fOT  one-dimensional 
silver  surfaces  with  o  =  X/lOO  and  t  *  X/5,  they  showed  that  a  peak  of  intensity  occurs  in  the 
backscatter  direction  for  p-polarisation  but  n<M  for  s-polarisation.  (In  p-polarisation,  the  electric 
vector  lies  in  the  plane  of  incidence,  whereas  for  s-polarisation  the  electric  vector  is 
perpendicular  to  Ae  plane  of  incidence:  p-polarisation  is  also  referred  to  as  TM  and  as  V, 
whereas  s-polarisation  is  also  called  TE  or  H.)  The  polarisation  dependence  suggests  that 
surface  plasmons  play  a  crucial  role  for  scattering  by  surfaces  of  these  parameters.  (Carefully 
controlled  experiments  to  verity  this  theory  have  still  to  be  made  but  Gu  et  al^  have  observed 
the  effect  for  a  (relatively  uncharacterised)  surface. 

Penurbation  theory  is  valid  only  for  relatively  smooth  surfaces  and  the  only  analytical  approach 
capable  of  dealing  with  rougher  surfaces,  o  •>  X,  is  physical  optics.  This  does  not  predict 
enhanced  backscattering  (see,  e.g.,  the  book  by  Beckinann  and  Spizzichino^).  In  experiments 
carried  out  in  1986,  inspired  by  reports  of  enhanced  backscattering  in  dense  volume  media  by 
three  groups''^,  Mendez  and  OTtonnelP-*  showed  experimentally  that  surfaces  which  show 
multiple  scattering  also  can  exhibit  enhanced  teckscattering* . 

The  work  of  M6ndez  and  ODonnell  has  proved  to  mark  a  watershed  in  progress  in  scattering 
from  randomly  rough  surfaces,  for  three  reasons.  Hrst,  particular  attention  was  paid  to 
fabricating  the  surfaces  with  reasonably  well-defined  statistical  properties;  the  light  scattering 
measurements  were  therefore  from  (more-or-less)  well  characterise  surfaces.  Second,  the 
scattering  cross-sections  were  remarkable  fw  that  fact  that  they  were  in  complete  disagreement 
with  existing  theoretical  predictions,  not  only  in  the  backscatter  direction  but  in  all  other 
directions  as  well.  Third,  the  impmtance  of  multiple  scattering  could  not  be  ignored.  Prior  to 
their  experiments,  enhanced  backscattering  due  to  scattering  fr^  high-sloped  surfaces  had  not 
been  pr^icted  by  a  single  theoretical  calculation  —  after  publication  of  the  results,  almost 
every  calculation  showed  the  effect! 

Work  in  rough  surface  scattering  at  Imperial  College  commenced  in  1985,  with  support  from 
the  UK  Science  and  Engineering  Research  Council.  During  the  first  year,  the  scattering 
equipment  was  designed  and  constructed  and  with  it  M6ndez  and  O'Donnell  obtained  their  first 


t  It  is  worth  pointing  out,  for  the  historical  record,  that  the  original  paper  by  Mendez  and 
O'Donnell  (ref  7)  was  rejected  by  Physical  Review  Letters  as  "not  of  interest  to  readers  of 
PRL",  despite  its  close  connection  to  the  work  on  dense  volume  media  which  was,  and  still  is, 
almost  exclusively  reported  in  PRL  and  PR. 
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enhanced  backscattering  results  in  mid- 1986.  On  completion  of  their  post-doctoral  terms, 
M6ndez  and  O'Donnell  took  up  permanent  positions  at  QCESE,  Mexico  and  Georgia  Tech, 
USA.  ARO  support  for  this  research  start^  in  mid-1987. 

Progress  during  period  July  ’87  to  July  '90 

One  dimensional  surfaces  -  experimental  results  for  conducting  surfaces 

The  experimental  results  reported  in  Refs[7]  and  [8]  were  for  two  dimensional,  isotropic, 
surfaces.  However,  the  prospects  for  theoretical  or  numerical  work  on  two  dimensional 
surfaces  were  not  good  and  since  a  major  aim  of  this  project  was  to  link  experiment  and  theory, 
it  was  decided  early  on  to  manufacture  and  make  measurements  on  one  dimensional  surfaces, 
i.e.  random  gratings.  The  surfaces  were  made  by  exposing  photoresist  to  laser  speckle 
patterns  produced  by  slit-like  apertures;  during  the  course  of  the  work,  replication  techniques 
were  developed  (see  [R2]).  As  with  previous  work,  all  surfaces  have  had  a  Gaussian 
probability  density  function  of  surface  height,  with  rms  roughness  o,  and  a  Gaussian 
correlation  function  of  l/e  width  equal  to  T.  Gt  is  possible  with  this  technique  to  make 
surfaces  with  a  probability  density  function  that  is  apptroximately  a  chi^  with  2N  degrees  of 
freedom,  where  N  is  the  number  of  independent  speckle  patterns  superposed  in  the  exposure  of 
the  photo-resist)  Another  advantage  of  one  dimensional  surfaces  is  that  they  can  be 
characterised  with  much  greater  confidence  than  two  dimensional  surfaces;  Ais  is  because  a 
sharper,  chisel-shaped  stylus  can  be  used  in  a  mechanical  surface  profiling  instrument  A 
Talystep  profiler  was  purchased  in  1988  for  this  project  using  University  funds  and  has 
significantly  unproved  the  reliability  of  surface  characterisation.  The  details  of  surface 
fabrication,  replication  and  measurement  are  given  in  Chapter  3  of  the  Appendix. 

The  emphasis  has  been  on  producing  a  few,  reliable  measurements  of  the  angular  scattering  by 
one  dimensional  rough  surfaces  at  up  to  four  wavelen^s  (X  =  0.63, 1.15,  3.39  and  10.6 
pm),  for  comparison  with  numerical  calculations.  Publications  [R3]  and  [R8]  give  "definitive" 
measurements  of  the  s-s  and  p-p  scattering  cross-sections.  In  [R8],  the  surface  has  an  rms 
roughness  a  =  1.22  ±  0.02pm  and  a  correlation  length  x  =  3.17  ±  0.07pm:  particular  attention 
should  be  given  to  Figures  3  and  4  in  [R8],  which  show  the  angular  scatter  for  p-  and  s- 
polarisation  at  three  angles  of  incidence  for  a  gold-coated  surface  at  the  four  wavelengths. 

These  results  are  available  on  floppy  disc  (PC  ot  MAC)  to  any  interested  person.  It  should  be 
noted  that  the  scattering  cross-sections  for  p-  and  s-polarisation  are  essentially  the  same  at 
shorter  wavelengths  (with  interesting  differences  in  the  side-lobes  of  the  enhanced 
backscattering  peak)  but  are  significantly  different  at  longer  wavelengths,  with  the  p- 
polarisation  showing  the  tendency  to  enhanced  backscattering.  This  is  consistent  with  the  fact 
that,  in  the  surface  plasmon  regime  (o«X),  wily  the  p-polarisation  couples  into  the  plasmons. 

We  have  compared  our  measurements  to  numerical  calculations.  These  have  been  done  more- 
or-less  independently  by  Nieto- Vesperinas^-*®,  MaradudinD.lZand  Maystre*^  for  the  kind  of 
surface  parameters  used  in  the  experiments.  Using  the  method  described  by  Nieto- 
Vesperinas^,  we  have  developed  programs  for  calculating  the  scanering  cross-section  for 
surfaces  that  are  (i)  perfea  conductors,  (ii)  non-absorbing  dielectrics  and  (iii)  arbitrary 
materials  (the  storage  requirements  and  computer  times  for  (ii)  twice  those  for  (i),  are  (iii) 
involves  another  doubling.)  In  [R3]  and  [R8]  we  compare  some  perfect  conductor  calculations 
with  experimental  results;  there  is  good  overall  agreement  at  near-normal  angles  of  incidence 
but  the  agreement  is  less  good  for  higher  angles  of  incidence  (>20°),  particulaiiy  at  higher 
scattering  angles  (>60°).  The  reason  for  this  discrepancy  is  still  not  clear.  It  should  stressed 
that  for  these  surface  parameters,  the  s-s  and  p-p  scattering  cross-sections  do  not  appear  to  be 
particularly  material  deptendent . 


One  dimensional  surfaces  -  experimental  results  for  dielectric  surfaces 

Dielectric  surfaces  scatter  less  than  reflecting  ones,  and  therefore  on  simple  grounds  would  not 
be  expected  to  show  multiple  scattering  effects  strongly.  We  have  been  able  to  make  dielectric 
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replicas  of  our  surfaces  and  therefore  cixnpare  the  scattering  cross-sections  of  identical  surfaces 
in  different  materials  ^old  and  silicone  rubber).  The  results  are  reported  in  [R3]  and  [R8]; 
perhaps  the  most  significant  feature  of  the  scattering  cross-sections  are  that  the  s-s  and  p-p 
sections  are  quite  different  and  very  recently  we  have  been  able  to  show  that  these  curves  can 
be  explain^  essentially  of  the  grounds  of  single  scattering!^.  As  with  the  gold-coated 
surfaces,  the  agreement  between  theory  and  experiment  is  fair,  but  there  is  sdll  significant 
disagreement  at  large  scattering  angles. 

If  the  rough  dielectric  surface  is  above  a  reflecting  plane,  then  analytical  and  numerical  [R8] 
studies  predict  a  sharp  enhanced  backscatter  peak.  With  careful  sample  preparation  (fOT  details, 
see  Appendix),  it  was  possible  to  make  a  very  thin  dielectric  sample,  rough  cm  one  side  and 
coated  with  gold  on  the  other.  The  experimental  results  are  repotted  in  [R8]  and  compared  to 
numerical  predictions,  with  good  agreement. 

Nieto-Vesperinas  and  Sanchez-Gil  predicted  numerically  an  unusual  result  for  the  light 
transmitted  by  a  rough  dielectric  film:  the  mean  angje  erf  refraction  roughly  equals  the  angle  of 
incidence  (rather  than  Snell's  law  being  obeyed).  This  can  also  be  explained  on  a  single 
scattering  argument  We  verified  their  prediction  experimentally  [Rrfj. 

Multiple  scattering  Kirchhoff  approximation 

The  Kirchhoff  approximation,  or  tangent  plane  approximation,  assumes  that  the  surface  is 
locally  flat  i-e.  that  the  local  radius  of  curvature  is  everywhere  much  larger  than  the 
wavelength.  It  is  almost  invariably  applied  in  a  single  scattering  way  but  there  is  no  reason 
why  it  cannot  be  extended  to  the  multiple  scattering  case.  This  was  done,  for  a  perfect 
conductor  in  [R5]  and  [R8],  and  for  the  more  general  case  in  Ref  14.  The  point  of  applying  the 
Kirchhoff  approximation  is  to  give  more  physical  insight  to  the  problem  than  the  "rigorous " 
calculation,  by  separating  single,  double,  triple....  scattering  terms. 

The  interesting  result  is  that  the  enhanced  backscatter  does  not  occur  for  the  single  scatter  case, 
i.e.  it  is  associated  with  multiple  scattering. 


Dynamics  of  volume  scattering 

A  small  amount  of  work  was  carried  out  on  the  temporal  behaviour  of  the  enhanced 
backscattering  peak  produced  by  dense  volume  media  [Rl].  In  this  case,  time  scale  probes 
path  length  of  scattering  in  the  media  (the  longer  the  path,  the  shorter  the  time  scale). 
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Photon  correlation  studies  of  light  scattered  by  e  dense  disordered  medium  st  end  around  the  backscatter  direction 
reveal  a  polariution  dependence  of  the  correlation  function.  This  is  related  to  the  polaruation  dependence  of  the 
mean  intensity  of  the  enhanced  backscatter  associated  with  coherent  cooperative  effects  in  such  media. 


There  has  been  considerable  interest  recently  in  the 
phenomenon  of  enhanced  backscattering  of  light  from 
disordered  random  media,*'*’  rough  surfaces,*®"’®  and 
atmospheric  turbulence.’*-”  The  phenomenon  itself 
is  of  general  interest  and  has  practical  consequences 
in,  e.g.,  remote  sensing,  but  much  of  the  motivation  for 
studies  in  this  area  comes  from  the  possibility  of  ob¬ 
serving  the  localization  of  light.” 

Figure  1  shows  the  basic  experimental  layout  for 
measurements  of  enhanced  backscatter.  Figure  2(a) 
shows  the  average  scattered  intensity  as  a  function  of 
angle  for  a  scatterer  consisting  of  a  10%  concentration 
of  0.46-Mm -diameter  latex  spheres  in  water  for  linearly 
polarized  incident  light  of  S14-nm  wavelength.  The 
average  intensity  of  the  copolarized  and  cross-polar¬ 
ized  scattered  light  is  different,  with  the  copolarized 
intensity  at  the  backscatter  typically  being  of  the  or¬ 
der  of  1.7  times  that  of  the  intensity  off  backscatter 
(e.g.,  at  20-30  mrad).  Figure  2(b)  shows  a  photograph 
taken  in  the  backscatter  direction  for  copolarized 
light.  It  is  now  established*~®  that  the  enhancement  is 
due  to  the  fact  that  the  scattered  wave  due  to  the  n 
scattering  events  associated  with  the  wave  vectors  k,, 
k| ,  k2, . . .  k/  (where  ki  and  k/  are  the  incident  and  final 
wave  vectors)  is  identical  to  the  reverse  path  in  the 
backscatter  (i.e.,  k/  =  -k,)  for  the  copolarized  compo¬ 
nent.  Since  forward  and  reverse  amplitudes  are  iden¬ 
tical  for  the  copolarized  case,  the  waves  interfere,  and 
an  enhancement  of  up  to  a  factor  of  2  is  possible  in  the 
backscatter  direction.  In  practice,  the  enhancement 
is  less  than  two  for  the  copolarized  case  because  of  the 
contribution  of  single  scattering  and  multiple  scatter¬ 
ing  involving  the  same  scatterer.  For  the  cross-polar¬ 
ized  component,  the  forward  and  reverse  paths  are 
partially  coherent,  with  the  degree  of  coherence  de¬ 
creasing  as  the  order  of  scattering  increases — in  effect, 
only  the  shortest  paths  give  any  enhancement.  The 
enhancement  factor  for  this  case  is  typically  wi.3  in 
the  above  experimental  arrangement. 

The  half-width  at  half-height  of  the  enhanced  back¬ 
scatter  peak  is  governed  by  diffraction  and  is  of  the 
order  of  X//*,  where  1*  is  the  effective  transport  mean 
free  path  length.  The  observed  backscatter  half-width 
of  the  order  of  3  mrad  in  Fig.  2(a)  implies  an  average 
transport  path  length  /*  «  170  ^m;  since  the  scattering 
mean  free  path  at  this  concentration  is  only  a  few 
micrometers,  it  is  clear  that  the  average  order  of  multi- 
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pie  scattering  is  very  high.  Because  of  diffraction 
effects,  all  path  len^hs  make  a  contribution  to  the 
copolarized  return  in  the  exact  backscatter  direction, 
whereas  only  the  shorter  paths  make  a  contribution  as 
the  angle  increases  away  from  backscatter,  since  the 
diffraction  cones  associated  with  shorter  paths  are  rel¬ 
atively  broader. 

The  dynamic  time  behavior  at  and  near  the  back¬ 
scatter  was  investigated  by  Maret  and  Wolf.’^  The 
correlation  function  is  strongly  nonexponential. 
Since  different  angles  of  observation  near  the  back¬ 
scatter  involve  different  total  transport  path  lengths, 
L,  each  of  which  has  a  different  time  scale,  the  tempo¬ 
ral  correlation  of  the  scattered  light  is  a  function  of 
the  angle.  The  field  correlation  is  given  by’^ 

(£(0)£*(t))  *  r  I{L)exp[-Lt/Wr)]dL, 

Jo 

where  HL)  is  the  intensity  contribution  of  transport 
paths  of  length  L  and  r  is  a  single  scattering  relaxation 
time.  Measurements  of  the  intensity  correlation, 

C(t)  =  (A/(0)A/(t)), 

are  given  in  Ref.  24  for  the  copolarized  case  and  are 
supported  by  the  above  reasoning. 

In  this  Letter  we  report  measurements  of  the  tem¬ 
poral  correlation  of  both  the  copolarized  and  cross- 
polarized  components  as  a  function  of  the  angle  near 
backscatter.  Correlation  times  are  extracted  as  a  de¬ 
scriptive  parameter  of  the  curves,  and  the  variation  of 
the  correlation  time  with  the  angle  close  to  backscatter 


Fig.  1.  Experimental  setup  for  measuring  the  intensity  of 
enhanced  backscattering  and  its  dynamic  properties.  BS. 
beam  splitter;  D,  diaphragm;  PMT,  photomultiplier  tube. 
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I  (a.u.) 


(a) 


(b) 

Fig.  2.  (a)  Measured  intensity  (in  arbitrary  unite)  as  a  func¬ 
tion  of  the  angle  around  backscatter.  Squares,  copolarized 
detection;  filled  diamonds,  cross-polarized  detection,  (b) 
Photograph  of  the  scattered  intensity  around  the  faackscat- 
ter  in  the  copolarized  case. 


turns  out  to  be  different  in  the  copolarized  and  cross- 
polarized  cases. 

In  our  experiments,  a  10%  concentration  of  0.46-pm- 
diameter  latex  spheres  in  water  was  used  with  an  ar¬ 
gon-ion  laser  at  X  =  514  nm.  The  sampling  pinhole 
diameter  of  approximately  30  prad  was  much  smaller 
than  the  speckle  size  (‘>'250  prad),  and  the  sampling 
time  of  2  psec  was  much  smaller  than  the  typical  corre¬ 
lation  time  of  100  psec.  A  polarizer  and  an  analyzer 
were  used  to  define  the  direction  of  linear  polarization 
of  the  illumination  and  detected  beams.  Standard 
photon  correlation  techniques  were  used  to  estimate 
the  intensity  autocorrelation  function  of  the  scattered 
light;  the  correlator  used  was  a  Langley-Fo.d  128- 
channel  instrument.  Typical  detected  photon  rates 
were  *0.5-1.0  per  sample  time,  and  a  typical  experi¬ 
ment  time  was  50  sec  (i.e.,  *'2.5  x  10''  samples). 

Figure  3  shows  a  typical  measured  correlation  func¬ 
tion  (for  the  copolarized  case  at  en  angle  S  -  3  mrad). 
The  general  form  of  the  curve  is  nonexponential,  in 
agreement  with  the  measurements  of  Maret  and 
Wolf.^^  The  observed  contrast  a*  (/)’  of  *>0.84  is  less 
than  the  Gaussian  speckle  value  of  unity  owing  to 
spatial  integration  and  dead-time  effects;  corrected 


for  these  effects,  the  normalized  contrast  is  approxi¬ 
mately  unity. 

Figure  4  summarizes  the  measured  correlation  data 
for  the  copolarized  and  cross-polarized  cases  at  angles 
of  0, 1, 3, 7,  and  9  mrad  from  backscatter;  the  vertical 
axis  is  the  logarithm  of  the  ratio  of  the  measured 
correlation  function  in  each  case  to  that  of  the  copolar¬ 
ized  case  at  a  large  angle  ( *40  mrad).  For  the  copolar¬ 
ized  case,  this  ratio  starts  at  approximately  unity  for  S 
~  0,  increases  to  a  maximum  at  0  *>  3  mrad,  and  then 
falls  to  unity  again  at  high  angles,  again  in  agreement 
with  previous  results.^*  However,  for  the  cross-polar¬ 
ized  case,  the  ratio  starts  at  unity  at  0  =  0  and  de¬ 
creases  smoothly  as  the  angle  increases. 

To  compare  the  differing  behavior  of  the  copolar¬ 
ized  and  cross-polarized  cases,  a  correlation  time,  de- 
flned  by 

,  r 

'  Jo  (iJ^) 

was  calculated  for  each  case  at  approximately  40  mea¬ 
surement  angles  from  -18  to  18  mrad.  The  results  are 


CIT) 


0  l,0E-04  2  0E-04  (S) 

Fig.  3.  Temporal  correlation  function  of  the  intensity  fluc¬ 
tuation  measured  at  0  =  3  mrad  for  the  copolarized  case,  with 
a  sampling  time  of  2  iisec  and  *0.6  detected  photons  per 
sample.  The  speckle  contrast  (r*/</)*  *  0.84. 


Lo9lO(T)/C|T)) 


Fig.  4.  Correlation  function  ratio  as  a  function  of  correla¬ 
tion  time  for  different  scattering  angles.  The  numbers  at 
the  right  are  the  different  measurement  angles  in  milliradi- 
ans  (numbers  without  parentheses  are  for  the  copolarized 
case;  numbers  with  parentheses  are  for  the  cross-polarized 
case).  C(71  is  the  temporal  correlation  for  0  »  40  mrad  in 
the  copolarized  case. 
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C.T.  (a.u.) 
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Fig.  5.  Measured  correlation  time  (C.T.)  in  arbitrary  units 
as  a  function  of  the  scattering  angle  for  the  copolarized 
(squares)  and  cross-polarized  (filled  diamonds)  cases. 

plotted  in  Fig.  5.  The  difference  in  behavior  of  the 
correlation  time  with  angle  for  the  two  cases  is  now 
clear,  with  the  copolarized  case  showing  a  double  peak 
and  the  cross-polarized  case  showing  a  gradual  de¬ 
crease  of  correlation  time  as  one  moves  from  exact 
backscatter. 

A  qualitative  explanation  of  the  observed  depen¬ 
dence  of  correlation  time  on  angle  and  polarization  is 
now  given,  and  it  is  based  on  the  fact  that  longer 
scattering  paths  give  rise  to  shorter  correlation 
times.^*  The  correlation  function  of  scattered  light 
consists  of  two  contributions,  one  from  the  incoherent 
intensity  and  the  other  from  the  coherent  intensity 
(i.e.,  that  responsible  for  the  intensity  enhancement  at 
backscatter).  The  incoherent  contribution  is  approx¬ 
imately  independent  of  angle  for  the  small  angles  of 
observation  around  backscatter  in  our  experiment, 
whereas  the  coherent  component  varies  with  the  angle 
and  polarization,  as  explained  below.  At  large  angles 
(i.e.,  a  few  tens  of  milliradians),  where  only  the  inco¬ 
herent  contribution  is  present,  the  time  scale  is  shorter 
for  the  cross-polarized  case  since  its  average  scattering 
path  length  is  longer. 

For  the  coherent  contribution,  two  effects  give  rise 
to  the  observed  behavior  of  the  correlation  time. 
First,  by  diffraction,  each  transport  path  length  con¬ 
tributes  an  enhancement  peak  whose  width  is  inverse¬ 
ly  proportional  to  the  path  length;  therefore,  observa¬ 
tions  at  8  =  0  include  coherent  contributions  from  all 
path  lengths,  whereas  observations  at  larger  angles 
include  coherent  contributions  only  from  shorter 
paths.  Second,  the  overall  magnitude  of  the  coherent 
contribution  decreases  away  from  the  backscatter,  and 
at  large  angles  the  time  scale  observed  is  that  of  the 
incoherent  component. 

For  the  copolarized  case,  it  has  been  established 
that  the  forward  and  the  reverse  scattering  paths  are 
always  coherent,  regardless  of  the  transport  mean  free 
path.^  All  transport  path  lengths  contribute  at  8  ~  0 
and  at  large  8,  giving  similar  correlation  times  in  these 
two  cases.  However,  as  9  increases  from  zero,  only 
shorter  transport  paths  contribute  to  the  coherent 
component  of  the  intensity,  and  thus  the  correlation 
time  begins  to  increase.  But  as  one  increases  8,  the 
relative  contribution  of  the  coherent  component  to  the 
total  intensity  decreases  and  ultimately,  at  large  an¬ 
gles,  the  correlation  time  must  decrease  again  to  that 


determined  by  the  incoherent  intensity.  The  maxi¬ 
mum  correlation  time  for  the  copolarized  case  is  ob¬ 
served  at  approximately  the  half-angle  of  the  enhance¬ 
ment  peak. 

For  the  cross-polarized  case,  only  the  shortest  for¬ 
ward  and  reverse  paths  are  coherent,*^  and  so  the 
dominant  effect  on  the  correlation  time  is  simply  the 
relative  contribution  of  the  coherent  component  to  the 
total  intensity,  which  gradually  decreases  from  a  max¬ 
imum  to  zero.  The  path-length  effect  is  negligibit,  as 
large  paths  never  contribute  to  the  coherent  comp.i- 
nent. 

In  conclusion,  we  have  shown  that  dynamic  light 
scattering  by  a  dense  disordered  medium  has  a  strong 
polarization  dependence  that  can  be  explained  quali¬ 
tatively  in  terms  of  the  proposed  mechanism  of  inten¬ 
sity  enhancement  shown  by  such  systems. 

This  research  was  supported  by  the  UK  Science  and 
Engineering  Research  Council  and  the  U.S.  Army  Eu¬ 
ropean  Research  Office. 
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A  method  of  replicating  randomly  rough  surfaces  fabricated  in  photoresist  has  been  developed,  enabling  compan> 
sons  of  the  scattered  light  to  be  made  between  identical  diffuaers  made  from  different  materials.  Experimental 
meaaurementa  on  one'dimenaional  metal  (gold)  and  dielectric  aurfacea  are  presented,  and  some  initial  compahaona 
with  numerical  calculationa  are  made. 


Recent  observations  by  Mendez  and  O’DonneU*-^  of 
enhanced  backscatter  and  strong  depolarization  from 
two-dimensional,  randomly  rough  surfaces  has  en¬ 
couraged  critical  discussion^  of  light-scattering 
mechanisms.  In  the  development  and  validation  of 
scattering  phenomena  the  provision  of  such  experi¬ 
mental  data  has  an  important  part  to  play.  In  this 
Letter  we  report  on  observations  of  light  scattering 
from  one-dimensional  randomly  rough  surfaces  at  a 
wavelength  of  0.633  >im,  comparing  a  metal-coated 
surface  with  a  dielectric  scatterer  having  identical  sur¬ 
face  characteristics;  some  preliminary  numerical  com¬ 
parisons  based  on  the  Kirchhoff-Helmholtz  integral 
equation'’  are  made  for  normal  incidence.  Few  experi¬ 
mental  results  have  been  reported  for  light  scattering 
from  dielectric  surfaces,^''^  least  of  all  any  that  com¬ 
pare  the  surface  scattering  properties  between  differ¬ 
ent  materials." 

The  requirement  for  producing  one-dimensional 
surfaces,  or  random  gratings,  arose  for  two  reasons; 
(1)  the  need  to  accommodate  one-dimenaional  scatter¬ 
ing  theories  and  (2)  the  computing  time  and  memory 
involved  in  running  numerical  codes  for  two-dimen¬ 
sional  problems. 

The  surfaces  are  produced  by  using  photoresist 
technology  similar  to  that  emploveil  i'y  Mendez  and 
O'Donnell’’^  but  using  a  thick  film  resist  (Shipley 
S1400-37)  to  give  a  film  thickness  of  ~12  iim.  The 
coated  substrates  are  exposed  to  eight  uncorrelated, 
one-dimensional  speckle  patterns  formed  by  focusing 
a  Gaussian  beam,  using  a  large-aperture  cylindrical 
lens  onto  a  ground-glass  diffuser.  The  resultant 
speckles  are  elongated  in  one  direction;  correlation 
lengths  of  the  order  of  micrometers  across  the  speckles 
and  of  millimeters  along  them  are  typically  achieved. 
Once  the  plates  are  rinsed  in  a  developer  whose  etch¬ 
ing  properties  are  linear  with  exposure  time,  the  one¬ 
dimensional,  randomly  rough  surface  is  completed. 

In  order  to  make  accurate  comparisons  between  the 
scattering  properties  of  different  materials  it  is  prefer¬ 
able  that  the  materials  have  identical  surface  charac¬ 
teristics.  By  using  the  etched  photoresist  plate  as  a 
master  it  is  possible  to  reproduce  the  surface  in  certain 
materials  by  a  method  of  replication.  Two  copies  of 
plate  #39  (rms  height  1.18  ±  0.13  fim,  1/e  correlation 


length  2.97  ±  0.05  /im)  were  formed  in  clear  silicone 
rubber  (Dow  Coming  Sylgard  182)  directly  from  the 
gold-coated  master  (the  gold  coating  was  applied  with 
evaporation  techniques  and  has  a  thickness  of  ~90 
nm).  The  first  copy  was  used  to  obtain  scattering 
measurements,  whereas  the  second  copy  was  used  to 
form  a  mold  from  which  an  epoxy  resin  replica  of  the 
original  was  cast  (Araldite  MY778  and  hardener 
HY956).  The  resin  copy  thus  forms  a  positive  replica 
of  the  original  surface.  In  order  to  determine  how 
successfully  the  master  had  been  replicated,  we  coated 
the  resin  copy  with  gold,  and  its  diffuse  scattering 
envelopes  were  measured  and  compared  with  those  of 
the  original  surface. 

The  scattering  equipment  and  geometry  used  are 

0*  incidence  pp 


Scaitotni  uif  k  /  depwi 


0”  incidence  o 


Sc«iienn(  antk  /defrees 

Fig.  1.  Comparuons  of  light  scattered  by  original  and  ep¬ 
oxy  resin  replica  of  plate  ^39  at  a  wavelength  of  0.633  ma. 
1/e  correlation  length  2.97  ±  0.05  /im,  rms  height  1.18  ±  0.13 
MXd;  both  surfaces  are  gold  coated. 
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Fig.  2.  AutocorreUtion  function  C(d)  «  (h(z)h(z  +  d))/ 
(hHz))  (where  (. . .)  represents  an  ensemble  average)  and 
probability-density  function  of  the  height  fluctuation  of  an 
epoxy  resin  replica  of  plate  #  39.  1/e  correlation  length  2.97 
±  0.05  pm,  rms  height  1.18  i  0.13  pm. 


essentially  identical  with  those  described  in  Ref.  2. 
For  all  scattering  experiments  the  incident  light  was 
linearly  polarized,  the  electric  vector  being  either  par¬ 
allel  (a  or  TE)  or  perpendicular  (p  or  TM)  to  the 
grooves.  The  results  are  shown  in  Fig.  1  for  normal 
incidence  only,  although  measurements  at  angles  as 
great  as  60  deg  have  been  carried  out.  The  scattereo- 
ligbt  envelopes  exhibit  the  phenomenon  of  enhanced 
backscatter  at  0  deg.  Each  measurement  was  normal¬ 
ized  assuming  perfect  conductivity  of  the  gold  coating, 
integrating  the  total  scattered  radiation  to  unity.  No 
depolarization,  either  sp  or  ps,  of  the  incident  radia¬ 
tion  was  observed,  as  expected  for  surface  roughness  in 
one  dimeimion  only.  The  scattering  properties  of  the 
replica  show  excellent  agreement  with  those  of  the 
original  surface.  The  slight  discrepancies  can  be  at¬ 
tributed  to  misalignment  of  the  scattering  and  detec¬ 
tion  planes  and  possible  fine-scale  resolution  limita¬ 
tions  («1  pm)  in  either  the  silicone  rubber  or  the 
epoxy  resin. 

Autocorrelation  and  height  probability-density 
functions  were  obtained  from  Talystep  measurements 
of  the  gold-coated  replica,  examples  of  which  are 
shown  in  Fig.  2.  Each  trace  consists  of  9000  data 
points  taken  every  0.2  pm  at  a  scan  speed  of  2.5  pm 
sec~'.  Six  uncorrelated  traces  displaced  across  the 
surface  (i.e.,  across  the  grooves)  were  averaged  to  ar¬ 
rive  at  the  quoted  characteristics.  The  results  in  Fig. 
2  show  good  agreement  with  Gaussian  distributions 
for  the  height  probability  and  autocorrelation  func¬ 
tions.  From  Fig.  1  the  exactness  of  the  agreement 
between  replica  and  original  indicates  that  Ae  initial 


copy  in  silicone  rubber  is  also  a  faithful  reproduction. 
Since  the  surface  statistics  are,  to  a  good  approxima¬ 
tion,  Gaussian,  the  fact  that  the  silicone  copy  is  a 
negative  of  the  original  should  not  affect  its  scattering 
properties. 

One  of  the  main  problems  involved  with  measuring 
the  light  scattered  frr  -n  the  surface  of  a  clear  dielectric 
medium  is  the  light  reflected  from  the  dielectric-air 
interface  at  the  back  of  the  sample.  The  unwanted 
transmitted  light  was .  bsorbed  with  a  neutral-density 
filter  of  density  4.0  optically  coupled  with  a  matching 
oil  (refractive  index  1.47)  to  the  Wk  face  of  the  repli¬ 
ca.  The  filter  was  angled  at  ~5°  to  the  vertical  to 
prevent  any  light  reflected  from  its  surface  from  enter¬ 
ing  the  detector.  The  diffuse  scattering  envelopes 
measured  for  the  dielectric  surface  are  shown  in  Fig.  3 
and  should  be  compared  with  the  relevant  graphs  in 
Figs.  1  and  4  of  the  perfect-conductor  measurements. 
Each  measurement  shows  the  diffuse  scattering  enve¬ 
lope  for  the  same  incident  power,  but  the  measure¬ 
ments  have  not  been  normalized,  and  the  intensity 
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Fig.  3.  Measurements  of  diffuse  scattering  envelopes  of  a 
dielectric  surface  illuminated  at  a  wavelength  of  0.633  pm. 
1/e  correlation  length  2.97  ±  0.05  pm,  rms  height  1.18  i  0.13 
pm.  Backscatter  occurs  to  the  right-hand  sides  of  the 
graphs. 
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Fig.  4.  Meaeuremente  of  diffuse  scattering  envelopes  of  a 
gold-coated  surface  illuminated  at  a  wavelength  of 0.633  ^m. 
1/e  correlation  length  2.97  ±  0.05  iim,  rms  height  1.18  0.13 
fim.  Backacatter  occurs  to  the  right-hand  sides  of  the 
graphs. 
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Fig.  5.  Numerical  calculations  of  diffuse  scattering  enve¬ 
lopes  for  a  dielectric  surface  of  refractive  indei  1.43  illumi¬ 
nated  at  normal  incidence  at  a  wavelength  of  0.633  am;  200 
realizations.  1/e  correlation  length  2.97  am,  rms  height  1.18 
am. 


scale  is  therefore  somewhat  arbitrary.  The  flrst  thing 
to  note  about  the  results  is  the  lack  of  any  signiHcant 
backacatter  enhancement  for  the  dielectric  surface,  in 
contrast  to  observations  for  the  gold-coated  surface. 
Both  ss  and  pp  data  have  similar  backscattered  inten¬ 
sities,  but  in  all  the  measurements  the  ss  data  show  a 
greater  reflectance  for  all  angles  than  the  pp  data.  A 
simple  explanation  may  be  found  by  considering  Fres¬ 
nel’s  reflection  formulas:  reflection  of  p-incident  ra¬ 
diation  from  a  plane  dielectric  surface  momentarily 
falls  to  zero  at  the  Brewster  angle,  whereas  reflection 
of  s-incident  radiation  rises  at  an  increasing  rate  with 
angle  of  incidence  and  is  always  greater  than  for  p- 
incident  radiation.  For  the  gold-coated  surface  (per¬ 
fect-conductor  case)  away  from  normal  incidence  the 
ss  and  pp  data  of  Fig.  4  show  only  slight  differences, 
whereas  the  corresponding  measurements  for  the  di¬ 
electric  case  are  quite  different  from  each  other.  On 
comparing  pp  data  for  the  two  cases  (metal  and  dielec¬ 
tric),  one  can  see  certain  similarities,  but  for  the  ss 
data  away  from  normal  incidence  one  measurement  is 
almost  a  reflection  of  the  other  about  0°. 

Figure  5  shows  the  result  of  a  Monte  Carlo  numeri¬ 
cal  calculation^  for  a  perfect  dielectric  of  refractive 
index  1.43,  rms  surface  height  1.18  pm,  and  correlation 
length  2.97  pm;  200  realizations,  each  with  300  sam¬ 
pling  points  along  a  length  of  25.31  pm  (40  wave¬ 
lengths),  were  used  in  the  calculation.  Comparison  of 
Fig.  5  with  the  upper  graph  in  Fig.  3  shows  broad 
agreement. 
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An  experimental  study  of  backacatter  enhancement  from  rough  surfacee  is  presented.  The  Stokes  parameters  of 
the  average  scattered  light  from  two-dimensional  rough  surfaces  show  the  presence  of  an  unpolarixed  component, 
which  lends  support  to  the  multiply  scattering  ray  model.  Experimental  data  from  one-dimensional  rough  surfaces 
are  compared  with  numerical  calculation. 


INTRODUCTION 

The  phenomenon  of  light  scattering  from  rough  surfaces 
(random  or  otherwise)  has  attracted  much  attention,  both 
experimentally  and  theoretically.  This  subject  is  of  particu¬ 
lar  importance  in  areas  that  involve  using  a  wave,  either 
acoustic  or  electromagnetic,  as  a  probe  to  observe  material 
and  surface  properties,  e.g.,  interpretation  of  radar  returns 
(from  the  surface  of  the  Earth  as  well  as  from  other  plane¬ 
tary  bodies)  and  noncontact  surface  characterisation. 

Recently  the  enhancement  of  scattered  light  intensity  in 
the  backacatter  direction  from  metallic  rough  surfaces  was 
reported  by  Mender  and  O’Donnell.''*  This  backacatter 
peak  from  other  random  rough  surfaces  had  been  previously 
observed;  a  sudden  increase  in  the  brightness  of  the  Moon 
when  it  approached  its  fullness  was  reported  os  far  back  os 
1924  by  Markov,*  while  Oetking^  reported  this  effect  to  be 
present  when  scattered  light  from  rocks,  as  well  as  certain 
reference  samples,  was  observed.  This  backacatter  peak, 
also  called  the  opposition  effect  in  the  literature,  has  been 
reported  by  several  other  authors*-*;  the  peaks  usually  have  a 
small  angular  width  (typically  2  to  3  deg)  and  result  when 
the  scattering  of  light  occurs  in  the  volume  as  well  as  on  the 
surface.  The  phenomenon  of  backacatter  enhancement  as 
observed  by  Mendez  and  O'Donnell  differed  from  what  the 
other  authors  reported  in  that  they  used  metallic,  high- 
sloped,  single-scale,  Gaussian,  random  rough  surfaces, 
whose  standard  deviation  of  surface  height  was  much  larger 
than  the  incident  wavelength.  The  key  property  was  that 
the  scattering  of  light  was  confined  to  the  surface  owing  to  its 
metallic  nature.  'Ihis  effect  was  normally  accompanied  by  a 
large  cross-polarized  component. 

The  high-sloped  nature  of  the  surface  meant  that  multiple 
scattering  was  a  significant  contribution  to  the  scattered 
light  The  methods  that  con  be  used  to  explain  this  phe¬ 
nomenon  analytically  ore  limited  because  of  the  restrictions 
imposed  on  the  available  scattering  theories.  Physical  op¬ 
tics*  cannot  be  used  since  it  accounts  only  for  single  scatter¬ 
ing  and  only  when  the  surface  structures  ore  much  larger 
than  the  incident  wavelength.  Analytical  multiple-scatter¬ 
ing  theories*-'*  utilizing  the  extended  boundary  condition 
ore  expressed  in  a  perturbation  series,  and,  owing  to  Uie 


difficulty  in  calculating  high-order  terms  and  the  slowness 
or  complete  lack  of  convergence,  they  have  been  limited  to 
the  case  of  low-sloped  surfaces.  The  full-wave  solution" 
may  account  for  the  enhanced  backacatter  peak  wbiU  sbigle 
scattering  only  is  used. 

The  technique  of  numerical  calculation  of  the  scattered 
light  has  been  available  for  some  time.'*-'*  This  method  is 
computationally  highly  intensive,  and  hence  the  application 
of  this  method  was  quite  limited  until  recently.  Perfectly 
conducting  surfaces  were  normally  considered,  and  the  in¬ 
tent  was  largely  to  establish  the  range  of  validity  of  the 
available  scattering  theories  such  as  the  physical-optics  so¬ 
lution'*'*  and  the  full-wave  theory.'*  The  technique  has 
also  been  used  to  calculate  the  scattered  light  from  high- 
sloped  surfaces,  and  enhanced  backacatter  peaks  were  ob¬ 
served  in  the  calculated  values'*;  such  effects  were  also  ob¬ 
served  when  real  metallic  surfaces  were  considered.'*  The 
major  flaw  of  this  procedure  is  that  it  gives  little  physical 
insight  into  the  scattering  process;  all  effects  such  as  shad¬ 
owing,  multiple  scattering,  and  light-surface  interaction  are 
mixed  inseparably  together. 

Mendez  and  O’Donnell  proposed  a  simple  model  involving 
multiple  scattering  of  rays  from  surface  facets;  it  was  an 
analog  of  that  used  in  the  volume-scattering  case.'*  Jake- 
man**  used  a  model  consisting  of  a  deep  random  phase 
screen  with  a  mirror  placed  just  behind  it  However,  al¬ 
though  these  models  explained  qualitatively  the  presence  of 
enhanced  backacatter  peaks,  they  leave  much  to  be  desired, 
e.g.,  they  do  not  take  into  account  the  polarization  of  the 
scattered  light  or  predict  the  detailed  shape  of  the  scatter 
envelope. 

The  research  reported  here  is  an  extension  of  that  present¬ 
ed  by  Mendez  and  O’Donnell.  The  normalization  method 
used  to  scale  the  experimental  data  will  be  discussed  briefly. 
The  polarization  behavior  of  the  scattered  light  from  two- 
dimensional  random  rough  surfaces  that  exhibit  enhanced 
backacatter  peaks  is  discussed.  Interpretation  of  Stokes 
parameters  leads  to  an  alternative  way  of  mapping  the  scat¬ 
tered  light;  instead  of  the  usual  copolarized  and  cross-polar¬ 
ized  intensity  one  can  plot  the  polarized  and  unpolarized 
components.  In  this  context,  unpolarized  means  that  there 
is  no  preferred  direction  of  polarization  over  the  solid  angle 
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of  measurement  (which  encompasses  many  speckles).  Fi- 
nally*  an  experimental  study  of  approximately  one-dimen¬ 
sional  random  rough  surfaces  is  presented,  and  the  experi¬ 
mental  data  are  compared  with  numerically  calculated  val¬ 
ues  where  possible. 

SCATTEROMETER  RESPONSE 

The  name  scatterometer  designates  the  equipment^'  that 
was  used  to  perform  experiments  involving  the  measure¬ 
ment  of  scattered  light  as  a  function  of  the  angle  of  incidence 
9,  and  the  scattering  angle  B. 

Figure  1  is  a  schematic  diagram  of  the  scatterometer 
viewed  from  above.  Two  laser  light  sources  were  available: 
He-Ne  (Spectra-Physics  105-1  laser,  wavelength  X  ■  0.633 
nm)  and  CO2  (Eklinburgh  Instruments  WL-4,  X  •  10.6  tim). 
The  spot  size  at  the  sample  was  approximately  10  mm  in 
diameter,  and  the  incident  beam  was  collimated.  The  sam¬ 
ple  mount  held  the  rough  surface  such  that  the  mean  surface 
normal  was  horizontal,  although  it  was  possible  to  tilt  the 
surface  normal  slightly  off  the  horizontal  plane.  The  sam¬ 
ple  mount  and  the  rotating  arm  were  movable,  both  having 
the  same  rotational  axis,  each  controlled  by  an  individual 
stepper  motor.  The  values  of  angles  were  such  that  was 
measured  clockwise  from  the  surface  normal  and  6  anticlock¬ 
wise  from  the  surface  normal,  ensuring  that  6  •  B,  in  the 
specular  direction.  The  detectors  used  were  a  Hamamatsu 
R647  photomultiplier  for  the  visible  light  and  a  Plessey 
PLT222  pyroelectric  detector  for  the  far-infrared  radiation. 
When  performing  experiments  in  the  far  infrared,  we  used  a 
chopper  in  conjunction  with  a  phase  lock-in  amplifier  (Stan¬ 
ford  Research  SR530)  to  eliminate  the  background  noise, 
the  incident  beam  being  chopped  at  typically  60  Hz.  A 
microscope  objective  of  5-mm  diameter  and  a  CdS  lens  of  1- 
cm  diameter  were  used  as  integrating  lenses  for  the  visible 
and  the  far  infrared,  respectively.  The  distance  between  the 
rotational  axis  and  the  integrating  lens  was  62  cm;  the  angu¬ 
lar  resolution  of  the  measurement  of  the  scattered  light  was 
thus  approximately  0.5°  for  the  visible  and  1.0°  for  the 
infrared. 

The  light  incident  upon  the  rough  surfaces  was  always 
coherent,  and  hence  a  speckle  pattern  was  generated.  The 
detector  response  is  proportional  to  the  spatial  integral  of 
speckles  in  the  solid  angle  of  the  integrating  lens,  i.e., 


PMT,  photomultipiier  tube. 
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Scattering  angle 

Fig.  2.  Scatter  envelope  from  a  MgO  surface,  s  incident  polariza¬ 
tion  (X  ■  0.633  «m).  Second  curves  from  bottom  denote  I,,.  curves 
immediately  below  these  indicate  Z,p,  and  topmost  curves  denote 
^1  is  a,  0*;  b.  — 30°,  and  c,  — 60*.  The  solid  curves  denote  the 
case  of  a  perfect  Lambertian  surface  with  a  perfect  scatterometer 
response. 


R^’-R  I  dnW(I?)J(n).  (1) 

J&Q 

where  AIJ  is  the  solid  angle  of  the  integrating  lens,  W(tl)  a 
weighting  function,  J?  is  the  constant  of  proportionality,  and 
J  is  the  radiant  power.  Here  we  make  an  assumption  that 
the  finite  spatial  average  is  equal  to  an  ensemble  average, 
i.e., 

(2) 

where  AO'  is  the  effective  integration  angle  and  does  not 
necessarily  stay  constant  for  all  6.  {J)  denotes  the  average 
of  Ute  radiant  intensity  for  that  particular  scattering  angle, 
for  an  ensemble  of  statistically  identical  but  independent 
rough  surfaces.  We  introduce  a  new  quantity,  the  mean 
normalized  differential  scattering  cross  section  (DSCS)  2, 
defined  as 
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where  is  the  incident  power.  £  is  related  to  the  well' 
known  bidirectional  reflection  function^^  if,)  by  a  simple 
expression: 

Z-  <A>cosfl.  (4) 

To  denote  the  polarization  property  of  the  incident  and 
scattered  light,  subscripts  are  appended  to  Z;  the  first  letter 
of  the  subscript  denotes  the  polarization  state  of  the  incident 
light,  and  the  second  the  polarization  state  of  the  detected 
light;  e.g.,  Z«p  U  the  DSCS  for  the  case  when  incident  light  is 
s  polarized  and  the  detected  light  p  polarized. 

Hence,  by  using  Eqs.  (2)  and  (3).  an  accurate  value  of  Z 
can  be  obtained  once  the  angular  dependence  of  AH'  is 
known,  and  this  is  usually  done  by  the  reference  sample 
method. Normally,  a  Lambertian  scatterer  is  sought  as  an 
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Fig.  3.  Scatter  envelope  from  MgO  lurfacc.p  incident  poiariution 
(X  ■  0.633  am).  Second  curves  from  bottom  denote  Zpp,  cui^s 
immediately  below  these  indicate  Xp,,  and  topmost  curves  denote 
Zi<K-  4,  isa,0*;b.-*30*,c. -60*.  The  solid  curves  denote  the  ease  of 
a  perfect  Lambertian  surface  with  a  perfect  scatterometar  retponae. 


Scattering  angle 

Fig.  4.  Plate  #313.  9,  »  0*,  X  »  0-633  ftm.  The  copolarized  and 
croas-poUrized  components  X„  and  Xgp  are  shown  in  a;  the  equiva¬ 
lent  Zpoi  and  ZuApo)  are  shown  in  b.  Solid  curves  denote  Z,,  and 
Zunpoi.  and  dotted  lines  denote  X,p  and  Zpo,. 

ideal  reference  sample.  If  the  directional  reflectance  (total 
scattered  power  divided  by  the  incident  power)  is  unity,  the 
mean  normalized  DSCS  for  a  Lambertian  diffuser  is  given 
by  a  simple  expression: 

lie,,  e)  =  --  cos  e.  (5) 

» 

The  problem  with  this  method  is  that  a  perfect  Lambertian 
diffuser  is  impossible  to  realize.  Historically,  a  freshly 
smoked  magnesium  oxide  (MgO)  surface  waa  used  as  an 
approximation  to  a  Lambertian  scatterer,  although  a  barium 
sulfate  (BaSO«)  surface  as  made  by  the  method  prescribed 
by  Eastman  Kodak  has  established  itself  as  a  standard.^* 
To  check  the  reaponae,  1  was  initially  obtained,  for  the  case 
of  a  MgO  surface,  assuming  that  Ad'  remains  constant  for  all 
e.  The  directional  reflectance  was  assumed  to  be  unity. 

Figures  2a,  2b,  and  2c  show  the  scattering  data  from  a 
MgO  surface  for  J,  «  O',  —30',  and  —60',  respectively.  The 
incident  light  waa  s  polarized,  of  wavelength  A  •  0.633  Mm. 
In  each  figure  the  copolarized  component  1„  and  the  croas- 
polarized  component  l,p  are  shown  together  with  the  turn  of 
both,  i.e.,  1m-  Figures  3a,  3b,  and  3c  show  similar  data  but 
with  the  incident  polarization  being  p  polarized.  The  fit 
with  the  cosine  line  is  fairly  good  for  both  incident  polariza¬ 
tions  except  for  the  presence  of  the  backacatter  peaks  (cf. 
Ref.  6  and  references  therein).  Since  our  scatterometar  gave 
an  approximately  coeine  reaponae  for  the  case  of  an  approxi¬ 
mately  Lambertian  scatterer,  it  was  decided  that  the  as¬ 
sumption  that  Alt'  remains  constant  for  all  f  was  a  reason¬ 
able  one  to  take.  This  approximation  is  implicit  from  this 
point  onward.  Another  important  approximation  used  waa 
that,  on  normalixation,  the  metallic  surfaces  were  all  as- 
•umed  to  be  perfect  conductors. 
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[.5.  Plate  #313,  » —10*,  X  »  0-633  Mtn.  Labels  as  in  Fig.  4. 
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Fif.6.  Plate  #313, », --20-,»- 0.633  »Di.  UbeU  u  in  Fif .  4. 
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STOKES  PARAMETERS  OF  SCATTERED  UGHT 

PUte  #313  was  a  gold-coated  metallic  two-dimensional 
Gauaaian  random  rough  surface.  It  was  made  by  the  meth¬ 
od  described  by  Gray,^  that  is,  by  multiply  exposing  a  pho- 
toresist-coated  substrate  to  laser  speckle  patterns.  After 
extensive  analysis  of  the  surface  profiles  obtained  from  a 
Talyaurf  profUometer  with  a  lufficiently  small  stylus  tip,  it 
was  found  that  plate  #313  had  a  standard  deviation  of 
surface  height  as  "  1.0  #  0.1  pm  and  a  1/e  correlation  length 
T  —  2.9  #  0.2  ant,  i.e.,  the  surface  structure  was  larger  than 
the  visible  ivavelength,  but  the  standard  deviation  of  the 
surface  slope  was  high,  approximately  0.5.  This  surface 
showed  the  enhanced  backscatter  peak  as  well  as  a  large 
crosa-polarixed  component.  In  fact,  the  scatter  envelope 
was  similar  to  that  obtained  from  plate  #  83  in  Ref.  2. 

Stokes  parameters  provide  a  complete  description  of  the 
polarization  state  of  the  scattered  light.  Ideally,  they 
should  have  been  obtained  for  all  t  and  for  all  angles  of 
incidence.  As  such,  they  were  calculated  for  normal  inci¬ 
dence  and  for  a  few  Fixed  scattering  angles  6  (one  near  the 
backscatter  direction).  The  incident  light  was  linearly  po¬ 
larized.  The  values  of  the  Stokes  parameters  showed  that 
for  those  fixed  angles,  the  scattered  light  was  composed  of 
two  components,  at  least  to  within  the  experimental  error 
(5%  of  the  first  Stokes  parameter):  one  that  is  li,.>.:;rly 
polarized  in  the  direction  of  the  incident  light  and  another 
that  is  unpolarized  (see  below).  For  scattering  of  monochro¬ 
matic  light  from  a  single  surface  realization,  light  scattered 
in  a  single  direction  will  be  completely  polarized.  However, 
what  we  are  measuring  is  the  scattered  light  integrated  over 
a  detector  solid  angle,  and  the  measurement  indicates  is  that 
there  is  a  component  of  scattered  light  that,  although  com- 
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Fig.  7.  Plate  #313,  -40*,  4 -0.633  sm.  Ubeli  as  in  Fig.  4. 
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Fig.  8.  Plate  #440,  s  incident  polariiation  (X  •  0.633  (.m).  for  »,  •  a,  0*;  b,  -10“;  c,  -20*;  d,  -40".  The  thicker  curves  denote  experimental 
X,,-,  the  thinner  curvee  ehow  the  numerically  calculated  valuea.  The  measured  croas-polarited  component,  denoted  by  the  line  almost 
coincident  with  the  x  axis,  is  shown  only  for  the  case  of  normal  tncideoce. 
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Fig.  9.  Plate  #  440,  p  incident  ^lariiation  (X  -  0.633  am),  for  -  a.  0*,  b,  -  lO";  c,  -20*,  d,  -40*.  Tha  thicker  curves  denote  eiperimenial 
ipp;  the  thinner  curvet  thow  the  numerically  calculated  valuea.  The  measured  cross-polarized  component,  denoted  by  the  line  almost 
coincident  with  the  x  azia,  it  shown  only  for  the  case  of  normal  incidence. 
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Fig.  10.  Piste  #440,5  incident  polarization  (A  »  10.6  ^m),  forS,  • 
a.  0*;  b.  -30*;  c,  “50*.  Labels  as  in  Fig.  9. 


pletely  polarized  in  a  single  direction,  does  not  have  a  pre¬ 
ferred  polarization  direction  averaged  over  the  detector  solid 
angle,  thus  giving  rise  to  a  measurement  that  makes  it  ap¬ 
pear  to  be  unpolarized. 

The  above  description  gives  rise  to  an  interesting  conjec¬ 
ture,  namely,  that  the  visible  scattered  light  averaged  over 
an  ensemble  of  rough  surfaces  such  as  plate  #313  is  com¬ 
posed  of  two  components  only,  that  is,  one  component  lin¬ 
early  polarized  in  the  direction  of  the  incident  polarization 
and  a  second  component  that  appears  to  be  unpolaiized. 
Keeping  this  in  mind,  it  is  thus  possible  to  describe  the 
scattered  light  not  in  terms  of  copolarized  and  cross-polar¬ 
ized  components,  say  1,,  and  Z,p,  but  in  terms  of  polarized 
and  unpolarized  components,  2poi  and  Xunpoii  respectively. 
The  relating  espressions  are 

2poi-Z«-2.p  <6a) 

and 


-  22,^.  (6b) 

Figure  4a  shows  the  plot  of  mean  normalized  DSCS  2w 
and  2.P.  while  Fig.  4b  ehowi  the  plot  of  mean  normalized 
DSCS  converted  in  the  manner  given  by  Eqa.  (6),  both  for 
the  case  of  plate  ^  313,  normal  incidence,  the  incident  wave¬ 
length’s  being  X  -  0.633  urn.  Figures  5.  6,  and  7  show  the 
similar  quantities  but  for  angles  of  incidence  of  -10*.  -20°, 
and  -40°,  respectively.  It  is  interesting  to  note  that  the 
enhanced  backscatter  peak  is  confined  mainly  to  the  unpo¬ 
larized  component. 

Multiply  scattering  rays,  according  to  the  model  proposed 
by  Mendez  and  O’Donnell,  are  the  cause  of  depolarization  as 
well  at  of  enhanced  backscatter  peaks,*  Le.,  they  will  not 
have  a  preferred  polarization  direction.  Singly  scattered 
rays  have  a  preferred  polarization  direction.  Hence  the 
unpolarized  component  will  contain  the  enhanced  backscat¬ 
ter  peak,  according  to  this  simple  model  This  is  the  case 
when  one  studies  Figs.  4-7. 
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Fig.  11.  Plate  #440,  p  incident  polarization  (X  *  10.6  am),  fort,  - 
a.0°;b.-30°;c,-60°.  Labels  at  in  Fig.  g. 
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Fig.  12.  PUtc  #440<X*  10.6  ^);  plot  of  the  aomnlimd  coherent 
component  power  for  $  end  p  incident  poUriution  obuined  from 
calculntion  and  eiperiment.  abown  aa  croaaea  and  open  circlea.  re- 
apectively.  The  physical-optica  solutions  for  a  and  p  polarizations 
and  a  one-dimenaional  rough  surface  arc  shown  dashed  and  dotted. 


SCATTERING  FROM  ON&DIMENSIONAL 
SURFACES 

Consider  a  surface  whose  surface  height  variation  depends 
on  one  Cartesian  coordinate,  i.e.,  a  one*dimensionai  surface 
2  «  h(x).  If  the  incident  electromagnetic  field  has  either  the 
electric  or  the  magnetic  field  lying  perpendicular  to  the 
incident  plane,  in  this  case  the  zr  plane,  there  is  no  cross- 
polarized  component,  and  the  scattering  problem  reduces  to 
a  scalar-wave  situation.  The  problem  still  remains  intracta¬ 
ble  analytically,  but  it  can  be  solved  exactly  by  numerical 
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calculation,  which  involves  averaging  the  calculated  scat¬ 
tered  intensity  over  an  ensemble  of  surface  realizations. 

Approximately  one -dimensional  random  surfaces  were 
made  by  etching  a  photoresist-coated  substrate  with  a 
speckle  pattern  whose  correlation  length  was  much  longer  in 
one  dimension  than  the  orthc^onal  one.  Two  surfaces,  both 
gold  coated,  were  considered:  plate  #440  *  1.2  ±  0.1 

Mm,T  *  2.9  ±  0.2  pm)  and  plate  #436  (at,  *  1.6  ±0.1  pm.r  » 
5.2  ±  0.3  pm).  Surface  parameters  were  obtained  by  using 
the  traces  measured  in  the  x  direction. 

Following  the  method  docribed  in  Refs.  15  and  16,  nu¬ 
merical  calculations  were  done  for  the  case  of  a  perfect  a>n- 
ductor  only.  To  ensure  the  accuracy  of  the  solution  ob¬ 
tained,  the  lengths  of  the  discrete  surface  (40X)  and  the 
sampling  distance  (0.13X)  were  chosen  such  that  the  fluctua¬ 
tion  of  the  normalized  total  scattered  energy  was  less  than 
3%. 

Figures  8a,  8b,  8c,  and  8d  show  for  6,  ®  0®,  —10®,  —20®, 
and  —40®,  respectively,  for  the  case  of  plate  #440,  the  inci¬ 
dent  wavelength  being  X  «  0.633  pm.  Figures  9a-9d  show 
2pp  for  the  same  angles  of  incidence.  Negligible  cross-polar- 
iz^  components  were  observed.  The  numerically  calculat¬ 
ed  values,  obtained  after  averaging  over  400  surface  realiza¬ 
tions,  are  shown  as  solid  curves.  Enhanced  backscatter 
peaks  are  observed,  both  in  the  experimental  data  and  in  the 
numerically  calculated  values.  Agreement  is  good  for  small 
6,  but  fails  when  $,  becomes  large.  The  reason  for  this  failure 
is  not  clear;  it  may  be  the  fact  that  the  finite  conductivity  of 
the  real  surface  has  not  been  taken  into  account,  that  the 
length  of  the  surface  in  the  computer  calculation  is  too  short, 
or  that  the  condition  that  states  that  the  total  scattered 
power  has  to  equal  the  incident  power  is  not  enough  to 
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Fig.  13.  Plate  #436  (X  ■  0.633  fim);  plots  of  Z„  (solid  curves)  and  Zpp  (dotted  curves)  for  9,  ■  a,  0®;  b,  -10®;  c,  -20®,  d,  —40®. 
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Fig.  14.  Plate  #436  (A  ■  0.633  #im);  plot  of  Z,,  (open  circles)  and 
(crosses)  for  B,  *  a.  0*;  b.  -30*;  and  c.  -50®. 


ensure  the  accuracy  of  the  solution,  as  is  known  from  other 
studies. 

Figures  10a,  10b,  and  10c  show  2„  for  6,  «  0®,  —30®,  and 
—50®  respectively,  for  the  same  plate  but  with  the  incident 
wavelength  now  being  X  «  10.6  pm.  Figures  lla>llc  show 
Ipp  for  the  same  angles  of  incidence.  We  believe  that  small 
cross-polarized  components  shown  in  Figs.  lOa  and  11a  are 
due  to  the  intrinsic  error  in  the  half-wave  plate  used  in  the 
scatterometer  and  the  imperfect  linear  polarization  state  of 
the  output  beam  from  the  COj  laser.  Numerically  calculat¬ 
ed  values  are  shown  for  comparison.  In  both  experimental 
data  and  numerically  calculated  values,  the  coherent  compo¬ 
nent  has  been  removed,  but  in  normalizing  the  experimental 
data  its  power  was  taken  into  account.  For  the  numerically 
calculated  values,  the  coherent  component  was  removed  by 
subtracting  from  the  total  intensity  the  modulus  square  of 
the  unscattered  field.^  The  agreements  are  good  for  all  9,. 
It  is  interesting  to  note  that  as  9,  increases,  the  total  power  of 
the  incoherent  component  seems  to  decrease  for  the  case  of  t 
incident  polarization,  while  for  the  p  incident  polarization  it 
seems  to  remain  constant.  This  becomes  more  evident 
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when  one  studies  Fig.  12,  which  plots  the  normalized  coher¬ 
ent  component  power  for  s  and  p  polarization,  experimental 
data,  and  calculated  values,  over  a  range  of  B,. 

Figures  13a,  13b,  13c,  and  13d  show  2„  and  Xpp  for  plate 
#436,  with  ■  0®,  -10®.  -20®,  and  -40®,  respectively,  for  X 
*  0.633  pm.  Figures  14a,  14b,  and  14c  show  Z„  and  Zpp  for 
the  same  plate,  with  8i  ■  0®,  -30®,  and  -50®,  respectively, 
for  X  •  10.6  pm.  Numerical  calculations  were  not  possible 
for  the  case  of  plate  #436,  as  the  required  array  sizes  were 
too  large  to  handle.  It  is  interesting  to  note  that,  unlike  for 
plate  #440,  there  seems  to  be  little  qualitative  difference  in 
the  scatter  envelopes  resulting  from  the  different  incident 
polarizations.  The  key  feature  to  note  is  the  small  or  nonex¬ 
istent  enhancement  in  the  backscatter  direction,  presum¬ 
ably  a  result  of  the  small  value  of  the  average  of  the  absolute 
slope. 

CONCLUSION 

Stokes  parameters  of  the  scatter  envelope,  which  shows  an 
enhanced  backscatter  peak,  have  been  studied;  it  was  found 
that  the  scattered  light  was  composed  of  polarized  and  unpo¬ 
larized  light  only  and  that  the  backscatter  peak  was  confined 
to  the  unpolarized  component  This  result  supports  the 
hypothesis  claiming  that  the  backscatter  peak  is  due  to  mul¬ 
tiple  scattering. 

An  experimental  study  of  scattering  of  light  from  approxi¬ 
mately  one-dimensional  surfaces  was  presented.  Compari¬ 
son  with  numerically  calculated  values  was  done  wherever 
possible.  There  are  disagreements  in  some  cases,  and  the 
reasons  are  unclear  at  this  stage  of  investigation.  Although 
numerical  calculation  gives  the  numbers  with  which  to  com¬ 
pare  the  experimental  data,  further  work  needs  be  done  if  we 
are  to  understand  the  physical  mechanism  behind  the  scat¬ 
tering  of  light  from  high-sloped  surfaces.  It  is  hoped  that 
the  experimental  results  presented  here  will  encourage  such 
work. 
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We  report  numerical  and  eiperimental  ofaaervationa  of  a  new  tranamiaaion  effect  in  rough  dielectric  interfaces. 


Recent  experiments  with  light  scattered  from  highly 
sloped  random  rough  surfaces'-^  and  the  development 
of  new  numerical  methods  of  study^'*  have  produced  a 
renewed  interest  in  scattering  effects  and  in  the  mech¬ 
anisms  that  produce  them,  in  both  metallic  and  dielec 
trie  8urfaces.^'“ 

We  report  the  numerical  prediction  and  experimen¬ 
tal  demonstration  of  a  new  transmission  effect  in  a 
semi-infinite  dielectric  material  when  light  is  incident 
from  a  semi-infinite  vacuum  upon  the  rough  surface 
that  separates  both  media. 

Because  of  the  limited  power  of  computer  calcula¬ 
tions  in  terms  of  both  speed  and  memory,  the  surface 
profiles  under  study  are  one  dimensional:  the  random 
height  depends  on  one  transverse  coordinate  only,  be¬ 
ing  constant  along  the  other  coordinate.  The  experi¬ 
ments  have  also  been  performed  on  specially  fabricat¬ 
ed  one-dimensional  surfaces  in  order  to  accommodate 
the  observations  to  the  available  theory. 

Numerical  Monte  Carlo  calculations  for  the  scat¬ 
tered  intensities  have  been  carried  out  from  samples 
with  rms  deviation  a,  normal  statistics,  and  a  Gaussian 
correlation  function,  with  profile  correlation  length  T. 
The  surfaces  are  generated  by  means  of  a  procedure 
identical  to  that  described  in  Refs.  3  and  7,  .namely, 
each  sample  consists  of  a  portion  of  a  sequence  of  ~10^ 
random  numbers  with  the  desired  normal  probability 
density  and  statistical  parameters  a  and  T.  The 
length  L  of  the  samples  is  typically  L  •  40  X  (A  being 
the  wavelength);  220  sampling  points  are  taken  for 
each  sample.  Both  the  field  and  its  normal  derivative 
are  calculated  on  the  surface  by  using  the  extinction 
theorem."  The  corresponding  equations  are  similar 
to  those  already  used  by  Maradudin  et  al.  in  Refs.  4 
and  10,  with  the  sole  exception  that  each  sample  is 
assumed  to  be  illuminated  by  a  plane  wave  instead  of  a 
beam  (i.e.,  no  tapering  is  made).  Then  those  bound¬ 
ary  values  of  the  fields  and  their  derivatives  on  the 
surface  are  introduced  into  the  far-zone  expression  of 
the  field  so  the  scattered  intensity  can  be  found.  Fi¬ 
nally  these  scattered  intensities  are  averaged  over  240 
samples.  At  a  given  angle  of  incidence  and  a  certain 
surface  record  of  length  L,  two  sets  of  scattered  inten¬ 
sities  are  obtained  by  considering  the  results  for  to  and 
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those  for  -So.  This  doubles  the  effective  number  of 
samples.  Osculations  were  done  on  a  CDC  Cyber 
180/855  computer.  Unitarity  of  the  normalized  re¬ 
flected  and  transmitted  intensities  at  the  dielectric 
interface,  as  well  as  the  convergence  of  the  results 
when  the  number  of  samples  is  increased,  is  taken  as  a 
criterion  of  numerical  consistency. 

The  angular  distribution  of  mean  intensity  in  the  far 
zone  of  light  transmitted  into  the  dielectric,  after  be¬ 
ing  scattered  at  the  rough  surface  separating  it  from 
vacuum,  is  plotted  in  Fig.  1 .  The  dielectric  constant  is 
SR(«)  =  1.991  [3(f)  *  0],  where  9}  and  3  denote  real  and 
imaginary  parts,  respectively,  and  the  surface  parame¬ 
ters  are  a  “  l.MX  and  T  »=  4.69X.  The  sum  of  the 
reflectance  H  and  the  transmittance  T  for  each  angle 
of  incidence  (0°,  20“,  40“,  and  60®,  respectively)  is 
shown  in  Fig.  1  (solid  curves,  p  polarization;  dashed 
curves,  s  polarization).  We  have  marked  the  direction 
of  refraction  that  would  correspond  to  a  plane  inter¬ 
face  (dotted  vertical  lines)  and  the  straight-through 
direction  (small  marks  at  the  left  of  the  vertical  lines  at 
the  tops  of  the  figures).  The  remarkable  effect  ob¬ 
served  is  the  deviation  from  the  refraction  direction  of 
the  mean  scattered  intensity  of  light  transmitted  into 
the  dielectric.  A  narrow  distribution  peaked  closer  to 
the  straight-through  direction  apjsears.  The  rough¬ 
ness  makes  the  angular  distribution  of  light  act  as 
though  there  were  no  different  refractive  index  on  the 
other  side.  Moreover,  this  peak  tends  to  grow  and  to 
become  narrower  as  the  angle  of  incidence  increases. 

We  believe  that  the  mechanism  that  produces  this 
effect  can  be  understood  within  the  diagrammatic  ap¬ 
proach  used  by  the  authors  of  Refs.  1,  2,  4,  9,  and  10. 
Because  of  the  slope  of  the  surface,  the  local  angle  of 
incidence  decreases  as  the  overall  angle  of  incidence 
increases.  On  the  other  hand,  since  the  material  is 
highly  transparent,  little  light  is  thrown  back  into  vac¬ 
uum  in  each  scattering  event,  so  double-scattering 
contributions  are  negligible  in  comparison  with  those 
of  single  scattering.  These  single-scattering  contribu¬ 
tions  tend  to  broaden  the  distribution  of  transmitted 
light,  but  the  dominant  angle  of  light  transmission  is 
observed  to  be  greater  than  it  would  be  if  there  were  no 
roughness  (Snell’s  law),  and  it  is  close  to  the  straight- 

C  1990  Optical  Society  of  America 


US2  OPTICS  LETTERS  /  Vol.  15,  No.  22  /  November  15, 1990 


through  direction:  thi»  is  due  to  the  aforementioned 
lower  local  angle  of  incidence.  In  the  case  of  noimegli- 
gible  multiple-scattering  contributions  (greater  re¬ 
fractive  index),  since  the  second  hit  on  ie  surface 
would  contain  almost  all  local  angles  of  incidence  the 
resulting  angular  distribution  should  be  expected  to 
spread  over  a  wide  solid  angle  failing  to  peak  clearly  in 
one  direction.  Calculations  made  by  increasing  i  have 
confirmed  this  point 

The  effect  reported  in  this  Letter  is  observed  nei¬ 
ther  on  surfaces  whose  radius  of  curvature  is  much 
larger  than  X  (this  is  the  general  criterion  of  validity  of 
the  Kirchhoff  approximation  for  perfect  conductors) 
nor  with  very  low  a  and  T  (say,  one  tenth  of  X).  How¬ 
ever,  as  this  effect  is  essentially  due  to  single  scatter¬ 
ing,  we  believe  that  a  calculation  based  on  the  Kirch¬ 
hoff  approximation  should  give  a  nearly  correct  solu¬ 
tion;  recent  preliminary  computations  confirm  this 
belief  (this  suggests  that  the  range  of  validity  of  the 
Kirchhoff  approximation  for  dielectric  surfaces  is 
much  broader  than  accepted  for  perfect  conductors). 

In  order  to  verify  the  fmdings  of  the  numeric^ 
transmission  calculations,  we  fabricated  a  dielectric 
diffuser  in  silicone  rubber  to  produce  a  slab  without  a 
tilt  angle.  The  diffuser  considered  here  is  a  dielectric 
replica  of  the  diffuser  investigated  in  reflection  by 
Sant  et  alfi 


Using  an  experimental  rig  essentially  identical  to 
that  described  by  O’Donnell  and  M6ndez,^  with  the 
exception  of  implementing  a  transmission  geometry, 
we  took  transmission  data  for  0°,  20°,  40°,  and  60° 
incidence  for  both  s  and  p  polarizations;  the  results  are 
shown  in  Fig.  2.  The  data  were  converted  by  Snell’s 
law  of  refraction  (for  the  plane  back  face  of  the  sam¬ 
ple)  to  represent  light  transmission  within  the  medi¬ 
um.  ’The  data  have  not  been  normalized  in  an  abso¬ 
lute  manner,  but  they  are  relatively  comparable. 

For  0°  incidence  the  transmission  peak  does  not 
appear  to  have  been  detected;  the  precise  reason  for 
this  is  not  known  but  is  thought  to  be  misalignment 
The  refractive  index  of  the  siHcone  rubber  at  X  =  633 
nm  has  been  more  accurately  determined  to  be  n 
1.411  since  it  was  reported  in  ^f.  6.  One  would  there¬ 
fore  expect  total  internal  reflection  to  occur  at  a  plane 
dielectric  interface  for  angles  of  incidence  beyond  arc- 
sin(l/1.411)  =*  ±45.1°:  this  is,  in  fact,  observed  for  40° 
and  60°  incidence,  where  the  detected  intensity 
abruptly  falls  to  zero  at  ~  —45°.  The  dashed  curve  on 
each  graph  of  Fig.  2  corresponds  to  the  refraction  angle 
at  which  the  transmitted  light  would  propagate  if 
there  were  no  front-surface  roughness. 

'The  agreement  between  the  experimental  and  nu¬ 
merically  calculated  results  is  excellent;  particular  at¬ 
tention  is  drawn  to  the  angular  position  of  peak  trans- 
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Fig.  1.  Numerical  calculations  of  transmitted  intensities  averaged  over  240  surface  samplea  of  length  L/X  -  40,  correlation 
length  r/X  «  4.69,  rms  deviation  a/X  ■  1.86,  and  dielectric  constant  e  “  1.991. 
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Fig.  2.  Mfc„3uremenU  of  transmission  envelopes  scattered  from  a  dielectric  surface  illuminated  at  a  wavelength  of  0.633  jim: 
1/e  correlation  length  2.97  ±  0.05  urn,  rms  height  1.18  ±  0.13  tan,  and  refractive  index  1.411. 


mission  and  angles  at  which  the  scattered  intensity 
falls  to  zero.  Although  for  60°  incidence  only  a  small 
portion  of  the  transmitted  light  is  detected,  agreement 
is  still  good. 
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Abstract.  Numerical  results  for  scattering  of  electromagnetic  waves  from  a 
rough  surface  are  presented  and  compared  with  experimental  results.  The 
method  uses  the  KirchhofT or  physical  optics  approximation  to  separate  the  single 
and  double  scatter  terms  in  the  total  scatter  pattern.  It  is  shown  that  enhanced 
backscatter  occurs  in  the  double  scatter  term  as  predicted  by  a  simple  ray  picture 
of  the  scattering  process. 

1.  Introduction 

The  observation  of  enhanced  backscattering  by  randomly  rough  surfaces  of  large 
root-mean-square  (r.m.s.)  slope  [1-4]  has  stimulated  a  re-examination  of  theoretical 
approaches  to  scattering  by  optically  rough  surfaces  [5-18].  It  is  now  possible,  at 
least  for  a  one-dimensional  surface,  to  solve  Helmholtz’  equation  numerically, 
realisation-by-realisation,  to  hnd  the  mean  scattered  intensity  and  its  statistical 
properties  (e  g.  variance).  This  is  a  powerful  method  for  predicting  the  angular 
distribution  of  the  scattered  intensity  and  has  been  exploited  by  Maradudin  (8-1 OJ. 
Maystre  [11],  Nieto-Vcsperinas  [12-14],  Thorsos  [17]  and  others.  Using  this 
method  one  can  predict  the  scattered  intensity  in  situations  where  controlled 
experiments  are  difficult  or  tedious. 

The  agreement  between  these  exact  calculations  and  experiment,  e.g.  as  reported 
in  [4],  is  good,  but  not  perfect.  The  problem  with  the  exact  method  is  that  it  give' 
little  physical  insight  into  the  scattering  process.  For  example,  the  exact  method  does 
not  answer  the  simple  question:  what  is  the  role  of  multiple  scattering  for  these  high- 
sloped  surfaces.^  In  contrast,  the  literature  on  light  scattering  by  dense  volume  media 
[19-21],  emphasizes  the  underlying  physical  processes  such  as  multiple  scattering. 

The  purpose  of  this  paper  is  to  use  an  approximate  theory  based  on  the  KirchhofT 
approximation  to  investigate  the  role  of  multiple  scattering  from  randomly  rough 
surfaces  of  large  r.m.s.  slope.  A  simple  picture  of  the  scattering  process  [1,2]  leads  to 
the  conclusion  that  the  backscatter  peak  is  caused  by  multiple  scatter  paths.  From 
figure  1  the  scattering  paths  ABC  and  CBA  give  scattered  waves  with  a  difference  in 
phase 

(k,  +  k2)*AC. 

This  term  is  zero  if 

kj^k,. 

which  is  exactly  the  case  for  backscatter.  Thus  the  backscatter  terms  add  coherently 
to  give  a  factor  of  two  larger  intensity  in  this  direction  than  in  other  directions  which 
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Figure  1 .  The  geometry  for  the  simple  ray  model  of  the  scattering  process  showing  a 
possible  path  ABC  and  its  time  reversed  partner  CBA. 

have  incoherently  added  terms.  The  width  of  the  backscatter  peak  can  also  be 
estimated  from  this  simple  model.  The  coherent  term  will  give  no  contribution  when 
the  scattered  waves  are  it  out  of  phase.  The  half  width  of  the  peak  is  then  found  to  be 

o _ L_ 

‘»~<|AC|>- 

To  calculate  the  held  scattered  by  a  surface^  first  the  boundary  conditions  on  the 
held  are  used  to  hnd  the  held  on  the  surface  and  then  the  scattered  (usually  far  held) 
distribution  of  this  held  is  found.  There  arc  two  categories  of  solution:  those  based  on 
exact  boundary  conditions  [5,6,8-16»22-28J  and  those  based  on  approximate 
boundary  conditions  [7,17,18,29-33]. 

When  exact  boundary  conditions  are  used,  the  solution  has  to  be  found  either 
numerically  as  in  [5, 6, 8-16],  or  approximately  using  an  appropriate  expansion  for 
the  scattered  held.  The  problem  then  is  that  it  is  generally  very  difficult  to  have  a 
good  physical  picture  of  what  the  equations  represent  since  the  expansions  of  the 
scattered  held  tend  to  be  rather  complicated.  This  means  that  the  methods  do  not 
give  much  insight  into  the  physical  mechanisms  behind  the  resulting  scattered 
intensity  distributions,  in  particular  whether  the  backscatter  intensity  is  due  to  single 
or  multiple  scattering  terms. 

The  second  group  of  methods,  in  particular  the  Kirchhoff  or  physical  optics 
methods,  can  however  give  a  clear  picture  of  the  physical  processes  involved  in  the 
scattering.  In  the  physical  optics  method  the  approximate  boundary'  conditions  at 
each  reflection  point  are  taken  to  give,  in  the  two-dimensional  case 

£,(x.ar)=(l+K)E,(x.a),  (1) 

~  »)• 


(2) 
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where  Ef^XyZ)  is  the  total  field  at  the  point  (*,«)  which  is  on  the  surface,  R  is  the 
planar  refiection  coefficient  at  that  point  which  depends  on  the  local  incidence  angle, 
£j(x,  a)  is  the  field  incident  at  that  point,  ki  is  the  incident  wave>vector  and  n  is  the 
outward  normal  to  the  surface  at  the  point  (x,  z).  This  simply  states  that  the  total  field 
on  the  surface  is  the  sum  of  the  incident  and  reflected  field  at  each  point. 

The  scattered  field  is  usually  taken  to  be  given  by  the  two-dimensional 
Helmholtz  integral  equation,  including  the  shadowing  functions 


dn 


(3) 


where  r  =  [(*  — is  the  zeroth  order  Hankel  function  of  the 
first  kind,  dr'  is  an  element  of  the  surface  and  and  represent  the 

incidence  and  scatter  shadow  functions  respectively. 


if  (x*,  a')  is  illuminated, 
if  (x',  x')  is  not  illuminated, 
if  (x',  x')  is  visible, 
if  (x',  x')  is  not  visible. 


It  is  important  to  note  that  these  are  geometrical  shadowing  functions  and  the  effects 
of  diffraction  are  not  taken  into  account.  The  use  of  these  straight  line  functions  is 
reasonable  if  the  distance  between  the  blocking  edge  and  the  shadowed  point  is  only  a 
few  wavelengths,  i.e,  if  the  shadowed  point  is  in  the  Fresnel  region  of  the  blocking 
point,  or  if  the  distance  between  the  source  point  and  the  blocking  point  is  of  the 
same  order.  This  situation  was  satisfied  in  the  cases  studied  with  the  distances  being 
of  the  order  of  the  correlation  length  of  the  surface  which  was  typically  3A. 

The  effects  of  the  shadowing  functions  were  discussed  by  Brown  [33]  who 
deduced  that  the  incidence  shadowing  function  was  valid  for  the  case  of  single  scatter 
whereas  the  scatter  shadow  function  showed  unphysical  behaviour.  Brown  sug¬ 
gested  that  the  scatter  shadow  function  causes  the  surface  field  to  depend  on  the 
position  of  the  observer  when  it  should  depend  only  on  the  incidence  angle  and  the 
form  of  the  surface  profile.  However,  S'(x',x')  does  not  change  the  surface  field,  it 
changes  the  observed  surface  field,  or  rather  it  changes  the  parts  of  the  surface  (and 
hence  the  surface  field  distribution)  which  are  visible.  This  is  not  unphysical.  The 
main  problem  with  ^fx',  x")  can  be  seen  from  figure  1 .  The  light  scattered  from  point 
A  in  the  direction  f*,  blocked  at  point  B  and  is  re-directed.  This  means  S(x\  x')  =  0 
and  the  contribution  from  this  light  is  neglected,  i.e.  the  re-directed  light  is  ignored. 
Therefore  S(x',  x')  in  equation  (3)  is  valid  for  the  case  of  single  scatter. 

Hence  investigations  of  the  validity  of  equation  (3)  are  therefore  investigations  of 
the  validity  of  two  separate  approximations. 


(i)  that  the  total  surface  field  and  its  normal  derivative  at  each  point  are  given  by 
equations  (1)  and  (2)  at  each  reflection  and 

(ii)  that  the  total  scattered  field  can  be  approximated  by  the  single  scatter  term  of 
equation  (3). 


In  this  work  only  the  first  approximation  is  used;  the  second  is  not  required  as  the 
second  order,  double  scatter,  term  is  explicitly  calculated.  This  involves  using  (i)  at  a 
scatter  point,  propagating  to  another  point  on  the  surface  taking  shadowing  into 
account,  using  (i)  again  at  the  second  scatter  point  and  then  propagating  to  the  far 
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zone  to  give  the  scattered  field.  This  term  ts  calculated  for  every  pair  of  points  to  give 
the  total  double  scatter  contribution. 

Obviously  this  second ^order  term  suffers  from  the  same  problem  as  equation  (3), 
in  that  there  will  be  a  scatter  shadow  function  which  will  predict  that  light  is  blocked 
and  re-directed  but  then  that  light  will  be  ignored.  However  the  energy  carried  by 
this  light  is  much  smaller  than  that  carried  by  the  double  scatter  term  and  it  is 
reasonable  to  expect  that  any  further  contribution  from  this  light  will  be  very  small. 
The  third-order  or  triple-scatter  has  also  been  calculated  for  a  few  cases  and  is  found 
to  be  negligibly  small. 

The  surfaces  considered  in  this  work  have  a  Gaussian  probability  density  of 
height  and  a  Gaussian  correlation  function.  Such  surfaces  are  manufactured  in  our 
laboratory  [3, 4]  using  the  technique  of  exposing  a  plate  with  a  layer  of  photo-resist  to 
many  independent  speckle  patterns  and  then  coating  with  gold.  The  experimental 
results  used  are  normalised  to  unit  area  beneath  the  curves  i.e.  the  material 
dependence  of  the  total  scattering  cross-section  has  been  removed. 

There  remains  the  question  of  validity  of  the  solution  obtained.  Equations  (1) 
and  (2)  contain  the  planar  reflection  coefficients  so  the  solution  should  be  valid  for 
surfaces  which  are  locally  quite  flat.  This  requirement  is  usually  interpreted  that  the 
radius  of  curvature  of  the  surface  should  be  very  much  less  than  the  wavelength  of 
the  incident  illumination.  The  radius  of  curvature  is  approximately  the  inverse  of  the 
second  derivative  of  the  surface  height  distribution.  As  the  height  distribution  is 
gaussian,  so  is  the  second  derivative,  with  a  standard  deviation 


where  tr^  is  the  standard  deviation  of  the  surface  height  distribution  and  t  is  the  1  le 
correlation  length.  One  of  the  surfaces  for  which  data  was  available  [4]  was  surface 
440  with  0^%  1*2  pm  and  r%  2*9  pm,  which  gives  -ftO’S  pm~  ^  and  thus  a  radius  of 
cur\'ature  of  2  pm.  Considering  illumination  with  0*633  pm  (He-Ne)  should  give 
a  case  where  most  of  the  radii  of  curvature  are  greater  than  the  wavelength  and  the 
approximation  (i)  is  reasonably  valid. 


2.  Theory 

The  single  scatter  term  is  given  by  equation  (3),  where  the  observation  point  has 
been  taken  to  be  in  the  far  zone.  Using  the  far  fields  of  the  Hankel  functions 


lim  H\l\kr)  = 


(4) 


and  equations  (1)  and  (2)  gives 


(5) 


—  (1  —  R)E,(,x',  z’)k  •  n  exp  ( —  ikj  ■  R)]  di’,  (6) 
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Figure  2.  The  geometry  used  for  the  Kirchhoff  approximation.  Note  that  the  incident  and 
scattered  angles  are  measured  in  opposite  senses. 


where  <p  is  A/?o“Jc/4,  Rq  is  the  vector  (», «)  and  R  is  the  vector  Also 

(see  figure  2). 

—  A  stn  —  k  cos 

k^ksinOx^kcos  Bt, 
n=  —  sin^K  +  cos^i, 


Then 


A  *  ft 

- 3*  tanp=s - =  m. 

cos  p  dx 


£■,(0)=— — 4^J* 

+  (I  —  R)k{m  sin  +  cos  0j)]£j(x',  z)  exp  ( —  ik  •  R)  dx'. 
The  incident  field  is  given  by  a  plane  wave 

Efx\  x')=.4jexp(iki  ■  R). 

Therefore,  following  Beckmann,  equation  (7)  reduces  to 

E,(e)  =  -  f— Y' Y  2  r 

•'  \nkTj  \  *  J  Lcos^  +  cosfl,  J 

X  J"  S{x',  »')5’(x',  z')RAi  exp  i(ki  —  k)  •  R  dx'  +  e(0), 


(7) 

(8) 


(9) 
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where  ^0)  is  an  edge  effect  term  which  is  much  smaller  than  the  first  term  inequation 
(9)  and  so  shall  be  neglected. 

Now  the  surface  is  discretised  into  N  segments  to  simplify  the  integration.  The 
phase  term  inside  the  integral  is  taken  to  be  constant  over  each  segment.  This  is 
simply  a  condition  on  the  number  of  segments  along  the  length  of  the  surface,  the 
length  being  small  enough  to  neglect  the  phase  term.  Taking  this  and  the  shadow 
terms  to  have  the  value  at  the  centre  of  each  segment  gives 


£,(0)=- 


2  exp(iy)  TH-cos(^-t-gi)~l 
nkr)  4  |_  cos^+cos^i  J 


N  Ci^ 

X  £/25(xj,a;)5'(x;,aj)expti(hi~k)-Ri)  A^dx\  (10) 

J-l  Ju-l)Ai 

with  the  surface  split  into  N  segments  of  equal  x-distance  Ax.  Assuming  a  small 
enough  value  of  surface  segment  length  Ax  gives 

fiax 

Aidx'=AiAx. 

J  (i-IIA* 

Considering  an  incident  plane  wave  of  unit  amplitude,  Ai=  1,  finally  gives 


E,(e)  =  ^ 


2  y^^fexpfiy)  ^ri+cos(g4-^i) 

Kkr  J  L  ^  |_  cos^+cos0j  _ 


X  J]/f5(x^,x})5'(x},x})exp[i(hi-h)*R^)Ax.  (11) 

This  is  the  usual  single  scatter  term  taken  to  be  the  solution  f^r  the  physical  optics 
approximation. 

Now  consider  the  double  scatter  contribution.  To  find  this  term  the  field 
obtained  by  scattering  the  incident  field  from  one  point  on  the  surface  to  another 
point  on  the  surface  must  be  found,  and  then  this  field  at  any  surface  point  due  to 
light  scattered  from  all  the  other  surface  points  is 


Jr  I  L  *^12  '^12  J 


-(t  -R^)^k^m,  sinS,  +  cos0i)HV*(fer,3)|-£i(2r,,*,)dX|,  (12) 

where  the  subscript  2  represents  the  second  point  and  subscript  I  the  first  point, 
r,2  =  [(*2  +  (*2“®  1  )*)'**• ‘he  normal  derivative  of  the  Hankel  function  is  given 

by 


Sn  L  'i2  'i!  J 


(13) 


and  the  shadow  function  5,]  is  defined  as 


if  (xj.aj)  is  visible  from  (x,,x,), 
if  (Xi,*,)  is  not  visible  from  (x,,x,). 
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The  field  scattered  from  all  these  second  points  will  then  be  the  double  scatter  term 
^xp(iy)  J 

+  i(l-R,)^^^^jexp(ikR,)dxj.  (14) 

Equations  (12)  and  (14)  together  form  the  double-scatter  term.  The  general  form  of 
this  term  can  be  seen  to  be  rather  complicated.  However,  if  p-polarised  (TM)  light 
incident  on  a  perfectly  conducting  rough  surface  is  considered,  for  both 

points  1  and  2.  Then  writing 

5^2“  *2)' 


E?K0)-- 


(  2  V''exp(i<p)  f  f  ^  r  *(*j-3c,)  *(22-a,)1 

-%kr)  .^J 


X  f/\*  *(fcr|  2)^(**>2  sin  6 — cos  0)  exp  [ifk^  —  h)  *  Rj]  dx2  dx , .  (15) 


For  the  discretisation  in  this  case  if  it  is  assumed  that  dx  is  small  enough  compared  to 
A  that  it  is  possible  to  simply  replace  each  integral  with  the  sum  of  the  values  at  the 
centre  point  of  each  segment  multiplied  by  Ax 


£121(0)  = 


'  2  Y'^exp(iy) 
itkr )  4i 


*(«|--2j) 

'0 


*(2|-2j)' 

'IJ  . 


X  W,'  \kT,j)k(m,  sin  0  -  cos  0)exp[i(ki-li)-RJAxAx.  (16) 
Similarily  for  s-polarisation  (TE).  £  =  —  1  at  both  points  and  t)ie  term  obtained  is 

I,  -  ^’] 


X  W,'*(*’'i;)*("'/S't>0i  +  co5  9j)  exp  [ilk,  — k)-  R,]  AiAx. 
To  obtain  this  result  the  following  approximation  was  used 

d 


(17) 


dn. 


J"  5(x, ,  2 1 )S ),k(m J  sin  + cos  d^EJ^Xj,  '(*r,p  dxj 


*  r  S(x„  2,)5^(mjSing)4-cose|)£|(2y,  Zf)  —  °  dx  (18) 
J  r  cn, 

i.e.  the  normal  derivative  of  the  shadow  function  Sj,  is  ignored.  The  effect  of  this 
approximation  on  the  final  result  is  not  known. 

Once  the  field  as  a  function  of  angle  is  known  the  power  per  unit  angle  at  the  far 
field  is  given  by 


J(0)  =  limr|£i(e)|2.  (19) 

r^oo 

For  a  perfect  conductor  /?  =  1  for  p-polarisation  and  -  1  for  s-polarisation,  therefore 
from  (11)  the  modulus  of  the  single  scatter  term  is  the  same  for  both  polarisations. 
The  double-scatter  term,  however,  is  different  for  the  two  cases,  ( 1 6)  for  p  and  (17) 
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for  s  polarisation.  The  sum  of  the  single  and  double  contributions,  the  'total'  power 
is  obtained  from 


7,(0)=limr|£'.(e)  +  Ci^'(e)l*. 

r-*** 


(20) 


In  this  case  the  sign  of  the  single-scatter  term  is  very  important  and  leads  to  a 
difference  in  the  total  scattered  intensity  distribution. 

A  test  of  the  reliability  of  any  method  for  calculating  the  scattered  field  is  the 
unitarity  or  the  power  scattered  divided  by  the  incident  power.  This  is  given  by 


P.(^)dy 

NdxcoiBi 


(21) 


3.  Results 

The  equations  for  the  single  and  double-scatter  terms  for  a  perfectly  conducting 
rough  surface  with  the  parameters  given  in  the  introduction  were  programmed  into  a 
computer.  To  reduce  speckle  noise  on  the  resulting  solutions  the  angular  distri¬ 
bution  of  intensity  was  averaged  over  many  (typically  1000)  different  realisations  of 
the  surface  profile  all  with  the  same  statistical  properties.  The  relative  amount  of 
energy,  the  unitarity,  in  the  single,  double  and  single  +  double  terms  is  shown  in  the 
table  for  both  polarisations.  From  this  cable  it  can  be  seen  that,  p^erhaps  surprisingl> . 
the  double-scatter  term  decreases  with  increasing  angle  of  incidence,  the  single 
scatter  increases  and  that  virtually  all  of  the  scattered  energy  is  contained  in  these 
terms.  It  also  appears  that  the  rate  of  decrease  of  the  amount  of  energy  in  thesccond- 
order  term  increases  with  increasing  angle,  for  s-polarisation. 

The  angular  distribution  of  the  scattered  light  is  shown  in  figures  3  and  4  and  the 
results  are  compared  with  experimental  cur\'es  in  figures  5  and  6.  Experimental 
results  are  represented  by  the  circles.  There  are  a  number  of  features  to  note  about 
these  results.  First  the  double-scatter  term  has  negligible  values  at  high  scatter 
angles  (higher  than  ±  60®)  so  this  is  the  region  where  the  single  scatter  dominates,  i.e. 
the  usual  single  scatter  physical  optics  approximation  value  is  closer  to  the  actual 
scattered  field  at  these  higher  angles. 


s-polartsation  (TE) 

Incidence 

Single 

Double 

Tot»l 

0“ 

0841 

0139 

0977 

10° 

0852 

0142 

0990 

20' 

0855 

ai35 

0987 

40° 

0886 

0088 

1  007 

p-polarisation  (TM) 

Incidence  Single  Double  Total 


0° 

0832 

0-156 

0-985 

10° 

0846 

0141 

0983 

20” 

0856 

0123 

0975 

40' 

a890 

0098 

1  014 
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Kifla/Artgto 


•e»Bw 


•can*r 


tcaiM'  angi* 


c  3.  Th«  single  (+  +  +),  double  (•••),  and  single  plus  double  ( — )  contributions  from 
surface  number  440  with  s-polarisation  incident.  The  backscattcr  angles  are  marked  by 
the  dashed  lines  to  the  right  of  the  graphs.  The  number  of  frames  averaged  over  are 
(fl)  2000  at  0%  (b)  1000  at  10'.  (r)  2000  at  20^  and  id)  1000  at  40^ 
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igure  S. 


Conin^ri>  i  of  single  plus  double  < — )  and  experimental  (OOO)  curves  for 
urface  number  440  with  sopolarisation  incident. 


Figure  6.  As  figure  5  but  with  p-polarisation  incident. 
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Secondly  there  is  no  narrow  peak  in  the  single  scatter  term  but  there  is  in  the 
double>scatter  one.  This  peak  occurs  at  the  backscatter  angle  and  there  are  signs  of 
side  lobes  on  either  side  of  this  peak .  Taking  the  intensity  at  the  maximum  of  the  peak 
and  dividing  by  the  intensity  at  the  same  angle  of  the  curve  obtained  by  simply 
drawing  a  smooth  line  joining  the  parts  of  the  curve  away  from  the  peak  gives  a  factor 
of  roughly  two.  This  is  the  factor  predicted  by  the  simple  picture  of  the  scattering 
process  involving  the  coherent  addition  of  waves  in  the  backscatter  direction.  It  is 
important  to  note  that  the  factor  of  two  enhancement  is  in  the  double  scatter  term  not 
the  total  term.  Again  this  is  predicted  in  the  simple  picture  which  requires  multiple 
scattering  to  have  coherent  interference  for  the  retro*reflected  light.  Obviously  this 
means  that  the  enhancement  factor  in  the  total  scatter  curve  will  be  less  than  two  in 
this  roughness  region  (t/tr^ 2%;  2*4)  as  the  single  scatter  is  not  zero  at  this  angle.  It  may 
be  that  the  factor  in  the  total  scatter  term  can  be  used  to  somehow  give  a  measure  of 
the  ‘roughness*  of  a  surface  since  as  the  roughness  is  increased  there  will  be  less  single 
scatter  and  more  multiple  scatter  giving  a  larger  enhancement. 

Another  point  which  is  apparent  is  that  since  the  amount  of  double  scatter 
decreases  and  the  double  scatter  curve  has  smaller  values  at  the  backscatter  angle  as 
the  incidence  angle  is  increased,  the  size  of  the  peak  and  hence  the  enhancement 
factor  decreases  with  increasing  angles  of  incidence. 


4.  Future  work 

Work  is  progressing  on  extending  the  method  to  the  triple-scatter  or  third-order 
term  and  also  to  the  double-scatter  term  from  dielectric  surfaces.  The  second-order 
dielectric  term  is  much  more  complicated  than  the  perfect  conductor  since  light  can 
travel  through  the  material  and  then  scatter  away  from  the  surface,  i.e.  light  can 
travel  between  grooves  which  would  be  hidden  from  each  other  considering  the 
usual  shadowing  term. 
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CHAPTER  15 


The  Opposition  Effect  in  Volume  and  Surface 
Scattering 
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1.  Introduction 

In  recent  years,  there  has  been  renewed  interest  in  the  so-called  opposi¬ 
tion  effect.  This  term  refers  to  the  enhanced  intensity  that  is  observed 
exactly  in  the  backscatter  or  retro-reflection  direction  when  certain  volume 
and  surface  scatterers  are  illuminated  with  parallel  light;  it  also  is  referred 
to  sometimes  as  enhanced  backscattering.  There  are  many  media  that 
exhibit  the  opposition  effect,  such  as  paper,  white  paint,  and  even  the 
surface  of  the  moon.  (The  moon  really  is  brighter  when  it  is  full  [1.2].) 
Materials  used  for  diffuse  reflectance  standards,  such  as  barium  sulphate, 
also  show  the  opposition  effect.  A  similar  phenomenon  can  be  seen  from 
scattering  by  clouds  (the  "glory”  or  "spectre  of  the  Brocken"  (3)).  al¬ 
though  the  mechanism  in  this  case  is  not  the  same  as  that  for  the  ex¬ 
amples  discussed  in  this  chapter.  Figure  I  shows  the  effect  for  the  case  of 
scattering  of  a  linearly  polarised  beam  by  a  dense  suspension  (about  10% 
by  volume)  of  latex  spheres  in  water.  The  upper  half  shows  a  photograph 
taken  in  the  backscatter  direction  and  the  lower  half  shows  the  co-  and 
cross-polarised  intensity  as  a  function  of  scattering  angle.  The  peak  inten¬ 
sity  in  the  backscatter  direction  in  this  case  is  up  to  twice  that  of  the  local 
background 

Observation  of  the  effect  is  well-documented  (4],  but  it  is  only  recently 
that  the  mechanisms  have  begun  to  be  understood.  In  this  chapter,  we 
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(b) 

Flo  I  Enhanced  backscauering  (opposition  effect)  by  a  dense  suspens.o?  ot  latex  spheres 
in  waier  (Fig,  2  of  D  N  Qu  and  J.  C.  Dainty.  Opr  Uu  13,  IO66-IO08.  19>^)  The  upper 
photograph  is  taken  in  the  backscatter  direction  (co-polansed  component)  and  the  lower  half 
show's  the  CO-  and  cross-polarised  intensities  as  a  function  of  the  scattering  angle 
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concentrate  on  one  particular  mechanism,  the  coherent  co-operative  effect 
of  multiple  scattering.  This  probably  is  the  dominant  mechanism  for  dense 
scatterers  like  white  paint,  but  in  other  cases,  there  are  geometrical  optics 
effects,  such  as  focussing,  that  also  can  give  enhancement.  To  give  some 
structure  to  the  chapter.  I  shall  discuss  the  cases  of  dense  volume  media, 
randomly  rough  surfaces,  and  double  passage  through  a  random  screen  in 
separate  sections;  the  key  unifying  feature  in  all  three  cases  is  the  occui- 
rence  of  multiple  scattering. 


2.  Scattering  by  Dense  Volume  Media 

The  scattering  by  many  surfaces  such  as  paper  and  paint  is,  in  fact, 
predominantly  a  volume  effect.  To  understand  the  scattering  process, 
experiments  often  are  carried  out  using  controlled  concentrations  of  mono¬ 
sized  latex  spheres  suspended  in  water  (S).  Although  the  opposition  effect 
has  been  documented  since  early  this  century,  in  recent  times  the  experi¬ 
ments  of  three  groups  [6-8]  have  resulted  in  a  renaissance  of  the  subject. 
(See  [9]  and  [iO]  as  a  starting  point  for  the  literature.)  The  reason  for  this  is 
not  so  much  the  importance  of  the  opposition  effect  itself,  but  that  the 
cause  of  it — multiple  scattering —  has  other  fascinating  effects — in  par¬ 
ticular,  the  possibility  ui  the  localisation  of  photons  in  a  disordered 
medium. 

The  basic  mechanism  for  enhanced  backscattering  is  illustrated  in  Fig.  2. 
If  a  linearly  polarised  plane  wave  is  incident  on  a  random  medium,  then  a 
ray  undergoing  multiple  scattering  with  particles  1.2.3 ...  A  has  precisely 
the  same  complex  amplitude  as  one  interacting  with  W. . .  3,2. 1;  the  two 
rays  are  coherent  and  thus  interfere  constructively,  giving  an  intensity  up 
to  twice  that  of  the  local  background,  where  there  are  no  such  coherent 
effects  A  more  rigorous  description  defines  each  scattering  event  by  a 


Fig  2.  Basic  multiple-scattering  explanaiion  of  the  opposition  effeci  forward  and  reverse 
rays  interfere  constructively  in  the  backscaiter  direaion. 
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matrix  equation,  in  which  an  incident  polarisation  vector  is  multiplied  by 
the  scattering  matrix  to  yield  the  polarisation  vector  of  the  scattered  light. 
Using  this  approach,  one  can  show  that  the  co-polarised  multiple-scattered 
light  is  enhanced  by  a  factor  of  two,  but  the  cross-polarised  forward  and 
reverse  paths  only  are  partially  coherent  and  give  an  enhancement  factor 
less  than  two:  in  fact,  as  the  degree  of  multiple  scattering  increases,  the 
enhancement  factor  tends  to  unity  for  the  cross-polarised  case. 

The  width  of  the  enhancement  peak  is  related  to  the  mean  transport 
length  of  photons  in  the  medium.  Long  paths  give  a  contribution  to  the 
enhancement  at  very  small  angles,  whereas  shorter  paths  contribute  over  a 
wider  range  of  angles.  Thus,  the  shape  of  the  enhancement  peak  gives 
information  on  the  length  of  light  paths  in  the  medium.  The  path  lengths 
also  are  contained  in  the  polarisation  behaviour,  since,  for  the  cross- 
polarised  case,  the  degree  of  correlation  between  the  forward  and  reverse 
paths  decreases  with  increased  order  of  scattering.  The  path  lengths  also 
can  be  probed  directly  using  femtosecond  pulses  and  observing  their 
broadening  after  scattering. 

The  energy  in  the  backscatter  peak  is  very  small  indeed  compared  to  the 
total  scattered  energy  in  all  other  directions,  so,  from  an  applied  point  of 
view,  one  could  reasonably  ask;  Why  is  there  all  this  fuss  (over  100 
publications  in  the  past  five  years,  symposia  and  workshops,  etc.)  over  such 
a  small  effect?  The  reason  probably  is  the  following:  The  opposition  effect 
is  the  result  of  scattering  mechanisms  that  hitherto  have  been  ignored 
almost  completely  and  that  could  have  profound  influence  on  our  under¬ 
standing  of  the  interaction  of  classical  waves  with  disordered  matter. 
Electrons  can  be  localised  in  a  medium  if  there  is  a  certain  degree  of  dis¬ 
order  and.  in  particular,  if  the  mean  free  path  of  the  electron  is  less  than 
A/2ir,  where  A  is  the  electron  wavelength.  If  one  could  make  such  a 
medium  for  visible  electromagnetic  waves,  light  inside  would  be  trapped, 
and  from  the  outside,  such  a  medium  would  be  a  perfect  reflector.  It  is  still 
not  clear  whether,  in  fact,  such  a  medium  can  exist  for  electromagnetic 
waves,  but  the  first  step,  that  of  weak  localisation,  has  been  observed  in  the 
manifestation  of  the  opposition  effect. 


3,  Scattering  by  Randomly  Rough  Surfaces 

There  is  an  enormous  amount  of  literature  on  the  elastic  scattering  of 
electromagnetic  radiation  by  randomly  rough  surfaces.  Because  an  exact 
general  solution  of  Maxwell's  equations  cannot  be  found  for  this  problem, 
the  theoretical  treatments  have  tended  to  focus  on  approximate  solutions, 
mainly  using  perturbation  theory  or  the  Kirchhoff  boundary  condition.  For 
surfaces  that  are  rough  compared  to  the  wavelength,  virtually  all  theoreti¬ 
cal  work,  until  recently,  has  assumed  single  scattering.  Experimental  work 
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on  rough  surface  scattering  has  suffered  from  the  fact  that  the  detailed 
surface  statistics  frequently  are  unknown  and  thus,  the  measurements  are 
only  of  empirical  value;  measurements  made  on  poorly  defined  surfaces 
are  of  no  value  for  checking  the  validity  of  approximate  theories. 

In  1987,  E.  R.  Mendez  and  K.  A.  O'Donnell  carried  out  some  of  the 
first  optical  scattering  measurements  on  well-defined  surfaces  with  relative¬ 
ly  simple  statistical  properties  111,12).  These  surfaces  were  made  by  expos¬ 
ing  photoresist  to  laser-produced  speckle  patterns  whose  statistics  are  well 
understood  [13|;  the  surfaces  can  be  coated  with  metal,  typically  gold,  or 
replicated  to  form  dielectric  or  metallised  copies  [  14).  Using  this  technique, 
surfaces  with  a  Gaussian  probability  density  of  surface  height  and  a  Gaus¬ 
sian  auto-correlation  function  can  be  produced  so  that,  for  the  first  time, 
surface  scattering  experiments  directly  relevant  to  theory  can  be 
performed. 

Some  of  the  first  experiments  conducted  by  Mendez  and  O'Donnell  were 
on  surfaces  of  root-mean-square  (rms)  surface  height  equal  to  approx¬ 
imately  two  wavelengths  and  correlation  length  five  wavelengths,  i  e  . 
surfaces  with  fairly  large  rms  slope  that  should  exhibit  multiple  scattering 
Figure  3  shows  one  set  of  angular  scattering  curves  for  a  one-dimensional 
version  of  such  a  surface  (15)  Note  the  opposition-effect  peak  in  the 
backscatter  direction.  This  has  the  same  origin  precisely  as  the  peak 
produced  by  scattering  from  dense  random  media;  forward  and  reverse  ray 
paths  are  coherent  and  thus  interfere. 


Scaitenng  angle 

Fio  3  Angular  scaucring  of  lighr  from  a  onr-dimensional  gold<oa«cd  randomly  rough 
surface  of  roor-mcan-square  height  fluctuation  a*  I.2i0.lMni  and  Gaussian  correlation 
function  of  1/e  length  2.9  2  0  2Mtn  (Fig  8  of  (15|)-  The  angles  of  incidence  arc  (a)  O', 
(b)  -lO®,  (c)  -20®,  and  (d)  -40®,  and  the  incidertt  light  is  j-polansed  (TE)  The  circles 
are  experimental  measurements  and  the  jagged  line  the  result  of  an  exact  calculation  for  a 
perfect  conductor  based  on  direct  numerical  solution  of  the  Helmholt^  equation.  Note  the 
enhanced  backscatter  peak  (opposition  effect)  in  the  backscatter  half-plane  on  the  right  hand 
side  of  each  graph. 


r  C.  Dainty 


Scattering  angle 


Scattering  angle 


Scattering  angle 


Fig.  3.  {Conunued) 
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Prior  to  the  experimental  work  of  Mendez  and  O’Donnell,  the  opposi¬ 
tion  effect  was  unknown  in  true  surface  scattering.  (Here,  we  ignore  the 
penetration  depth  of  the  wave  in  a  metal,  which  is  on  the  order  of  80A  for 
visible  light  and  a  gold  surface.)  Their  work  has  stimulated  a  large  number 
of  groups  to  search  for  a  more  “rigorous”  explanation  of  the  enhanced 
backscatter  peak;  electromagnetic  theorists  are  not  very  fond  of  rays!  Un¬ 
fortunately,  no  analytical  solution  exists  and  thus  one  has  to  resort  to  — 
or  exploit — numerical  techniques.  These  give  quite  good  agreement  with 
the  measurements  [15-18],  but  otherwise  give  little  physical  insight  to  the 
nature  of  the  scattering.  Figure  3  shows  the  result  of  an  "exact" 
numerical  calculation  for  a  perfect  conductor  based  on  the  Helmholtz 
equation.  The  agreement  is  fair,  although  there  is  a  significant  discrepancy 
at  higher  angles.  (This  almost  certainly  is  not  due  to  the  fact  that  the 
calculations  are  for  a  perfect  conductor,  whereas  the  experiments  were 
done  on  gold-coated  surfaces.) 

An  alternative  approach  is  to  adapt  a  theory  based  on  approximate 
boundary  conditions,  such  as  the  Kirchhoff  or  physical  optics  theory,  to 
multiple  scattering.  Figure  4  show-,  the  result  of  such  a  calculation  carried 


Fig.  4.  The  single  (  +  +  +  ).  double  (*••)  and  single-i-double  ( — )  coniribunons  lo  ihc 
angular  disiribulion  of  scanered  light  for  s-polarisation  mcidcnl.  for  the  same 
surface  parameters  as  those  in  Fig  3.  calculated  using  a  double-scattering  version  of  the 
Kirchhoff  approximation. 
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Fic  S  Photographs  of  the  backscatiered  iniensii)  from  a  two-dimensional  randomly 
rough  surface  taken  (a)  with  no  analyser,  (b)  through  an  analyser  parallel  to  the  incident 
polarisation,  and  (c)  with  the  analyser  perpendicular  to  ihe  incident  polarisation  (courtesy  of 
E  R  Mendez  and  K.  A.  O’Donnell) 
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out  to  a  double-scattering  approximation  (including  shadowing)  for  a  per¬ 
fectly  conducting  surface  of  the  same  parameters  as  that  in  Fig.  3,  illumin¬ 
ated  by  s-polarised  (TE)  light  of  wavelength  0.63  /im  at  normal  incidence. 
The  enhanced  backscatter  peak  shows  only  in  the  double-scatter  term, 
providing  strong  support  of  the  ray  theory  that  emphasises  the  crucial  role 
of  multiple  scattering. 

One  important  difference  between  the  surface  and  volume  scattering 
cases  is  that,  in  the  surface  case,  only  low-order  multiple  scattering  occurs, 
mainly  just  second  order.  This  means  that  the  co-  and  cross-polarised 
contributions  to  the  enhanced  backscatter  peak  are  equally  large.  For 
two-dimensional  rough  surfaces,  there  are  strong  polarisation  effects. 
Figure  S  shows  photographs  of  the  backscattered  light  for  a  surface  illumi¬ 
nated  by  linearly  polarised  light  and  viewed  (a)  with  no  analyser,  (b)  with 
the  analyser  parallel  to  the  incident  electric  vector,  and  (c)  with  it  perpen¬ 
dicular  to  the  incident  electric  vector.  The  effect  is  very  striking  when 
observed:  as  the  analyser  is  rotated  through  90°,  the  scattering  pattern 
appears  to  rotate  4S°.  Again,  a  simple  ray  argument  seems  to  explain  the 
observations  [12]. 

As  with  the  case  of  scattering  by  dense  volume  media,  it  is  not  the  op¬ 
position  effect  itself  that  is  particularly  important —  although  it  is  a  fasci¬ 
nating  and  unusual  phenomenon — but,  more  significantly,  the  fact  that  it 
encourages  us  to  question  the  accepted  treatments  of  a  subject.  For  exam¬ 
ple,  one  hears  reference  to  Lambertian  scattering  surfaces,  but  it  seems 
certain  now  that  such  a  thing  simply  is  a  hgment  of  the  imagination  and  not 
realisable  by  any  real  surface,  because  of  multiple  scattering.  Likewise,  the 
angular  scattering  curves  of  Fig.  3  are  quite  remarkable,  even  ignoring  the 
enhanced  backscatter  peak;  for  example,  the  existence  of  a  large  compo¬ 
nent  of  light  scattered  in  the  backscatter  half-space  at  angles  greater  than 
the  incident  angle  is  unexpected  and  counterintuitive. 


4.  Double  Passage  through  a  Random  Screen 

There  is  an  extensive  amount  of  literature  on  the  double  passage  of  light 
through  turbulence.  (See  [19]  for  a  partial  list.)  By  analogy  with  the  above 
cases  of  surface  and  volume  scattering,  one  expects  to  observe  the  opposi¬ 
tion  effect,  and  this  has  been  reported  (20,21).  A  practical  consequence  of 
this  is  that  diffraction-limited  information  is  present  in  the  average  image 
of  a  deterministic  object  that  is  illuminated  and  viewed  through  the  same 
distorting  medium.  As  an  example  of  this,  we  consider  the  simplest  case, 
that  of  a  Michelson  or  two-slit  interferometer  illuminating  and  viewing  an 
object  through  a  random  (phase)  screen,  as  in  Fig.  6. 
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Again,  we  consider  a  simple  ray  approach  to  explain  why  difFraciion- 
limited  information  is  present  in  the  image  of  a  point  object.  Consider  light 
transmitted  through  slit  A  towards  the  object.  It  accumulates  a  random 
phase  0,  on  its  passage  through  the  random  screen  and  then  is  scattered  in 
all  directions  by  the  point  object.  That  portion  of  the  light  that  returns 
through  slit  A  has  a  random  phase  26|  and  does  not  contribute  to  the 
diffraction-limited  process.  However,  another  portion  of  the  light  travels  to 
slit  B,  accumulating  a  total  random  phase  of  6,  Now  consider  light 
that  is  emitted  from  slit  B.  This  also  will  accumulate  a  total  random  phase 
oi  0i  +  $2  by  the  time  it  arrives  at  slit  A  and  thus,  constructive  interference 
will  occur  between  the  forward  (A  to  B)  and  reverse  (B  to  A)  waves.  This 
constructive  interference  is  superimposed  on  a  background  level  caused  by 
the  other  paths  A  to  A  and  B  to  B. 

Figure  7  shows  a  computer  simulation  of  the  fringes  in  a  Michelson 
interferometer  viewing  a  point  source  for  (a)  no  random  screen,  (b)  double 
passage  through  different  random  screens  (2.000  realisations),  and  (c) 
double  passage  through  the  same  random  screen;  (d),  (e),  and  (f)  show  the 
Fourier  transforms  of  (a),  (b),  and  (c),  respectively.  Note  the  high  contrast 
fringes  formed  in  (c)  and  corresponding  side-bands  in  (Q.  of  magnitude 
one-half  that  of  the  no-random-screen  case. 

This  technique  can  be  generalised  to  non-redundant  and  filled  imaging 
pupils  (19, 22],  with  corresponding  reductions  in  fringe  visibility. 


5.  Discussion 

The  increased  understanding  of  the  opposition  effect  over  the  past  five 
years  has  had  an  impact  in  at  least  three  areas  of  physics.  First,  the 
importance  of  multiple  scattering  has  been  exposed:  this  has  a  profound 
importance  for  the  behaviour  of  light  in  dense  scattering  media,  and  may 
lead  to  the  phenomenon  of  photon  localisation  in  disordered  structures. 
Dense  random  media  are  a  new  type  of  scattering  that  may.  in  fact,  be  the 
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(b)  Double  passage  •  different  screen 


(e)  Fourier  Transform  of  (b) 


(c)  Double  passage  through  same  screen 

(0  Fourier  Transform  of  (c). 

Fig.  7  Computer  simulation  of  the  fringes  formed  m  a  doublc'passage  Michelson 
interferometer 


218 


).  C.  Dainty 


first  natural  random  structure  to  show  the  unusual  spectral  shifts  predicted 
by  E  Wolf  (23.241. 

Secondly,  the  opposition  effect  from  well-defined  surfaces  has  led  to 
renewed  theoretical  studies,  mainly  using  numerical  techniques.  Consider¬ 
ing  that  this  effect  for  true  surface  scattering  was  not  predicted  prior  to  its 
experimental  discovery  in  1987,  it  is  quite  amusing  to  observe  that  almost 
every  study  of  light  scattering  is  now  considered  incomplete  without  a 
discussion  of  it!  Thirdly,  it  is  becoming  clear  that  the  opposition  effect — or 
more  precisely,  multiple  scattering — has  applications  in  areas  other  than 
scattering  in  volumes  and  surfaces.  The  example  of  the  double  passage  of 
light  through  a  random  screen  is  one  illustration  of  this,  and  doubtless 
there  will  be  others  in  the  future. 
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Abstract.  The  enhanced  backscattering  of  light  from  a  random  surface  is  manifested  by  a 
well  defined  peak  in  the  retrc-refieciion  direction  in  the  angular  distribution  of  the 
intensity  of  the  incoherent  component  of  the  light  scattered  from  such  a  surface.  In  this 
paper  we  present  several  new  theoretical  and  experimental  results  bearing  on  the 
conditions  under  which  enhanced  backscattering  occurs,  and  the  way  in  which  this 
phenomenon  depends  on  the  nature  of  the  random  surface  roughness,  both  in  the  case 
that  the  random  surface  bounds  a  semi-infinite  scattering  medium  and  in  the  that  it 
bounds  a  film,  either  free-standing  or  on  a  reflecting  substrate.  In  addition,  we  present 
new  results  on  the  transmission  light  through  thin  metallic  films  bounded  by  random 
surfaces,  which  display  the  phenomenon  of  enhanced  transmission,  namely  a  well  defined 
peak  in  the  antispecular  direction  in  the  angular  distribution  of  the  intensity  of  the 
incoherent  component  of  the  light  transmitted  through  such  films. 


1.  Introduction 

Enhanced  backscattering  in  (he  scattering  of  light  from  a  random  surface  is  manifested 
by  a  well  defined  peak  in  the  retro-reflection  direction  in  the  angular  dependence  of  the 
intensity  of  the  incoherent  component  of  the  light  scattered  from  such  a  surface.  It  has 
now  been  studied  intensively  theoretically  [1-15]  and  experimentally  [16-20].  On  the 
basis  of  these  studies  it  is  now  believed  that  the  enhanced  backscattering  of  light  from 
a  moderately  rough,  reflecting,  random  surface  results  primarily  from  the  coherent 
interference  of  each  doubly  reflected  optica]  path  with  its  time-reversed  partner.  If  the 
scattering  surface  supports  surface  electromagnetic  waves,  enhanced  backscattering 
occurs  even  for  weak  corrugations  [1-3,6.15].  In  this  case  it  is  due  primarily  to  the 
coherent  interference  of  each  doubly  scattered  surface  electromagnetic  wave  path  with 
its  time-reversed  partner.  Such  enhanced  backscattering  phenomena  are  examples  of  a 
broader  class  of  multiple  scattering  phenomena  that  go  under  the  name  of  weak 
localization,  in  the  present  case  of  light  and  of  surface  electromagnetic  waves,  respectively. 
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In  this  paper  we  present  new  computer  simulation  and  experimental  results  that 
bear  on  the  enhanced  backscattering  of  light  from  random  one-dimensional  surfaces  on 
semi-infinite  substrates  and  thin  films.  The  emphasis  in  the  investigations  described  here 
is  on  elucidating  the  mechanisms  responsible  for  enhanced  backscattering.  and  the  way 
in  which  this  phenomenon  depends  on  the  nature  of  the  surface  roughness 

In  addition,  we  present  new  computer  simulation  and  experimental  results  concerning 
the  enhanced  transmission  of  light  through  thin  metal  hims,  one  or  both  of  whose 
surfaces  are  random  one-dimensional  surfaces.  Enhanced  transmission  is  a  well  defined 
peak  in  the  antispecular  direction  in  the  angular  distribution  of  the  intensity  of  the 
incoherent  component  of  the  light  transmitted  through  such  films.  Of  interest  is  the  way 
in  which  this  phenomenon  depends  on  the  nature  of  the  surface  roughness. 

Thus,  we  will  be  concerned  with  the  interaction  of  light  with  surfaces  defined  by  the 
equation  X3  =  ((x,),  so  that  their  generators  are  parallel  to  the  Xj-axis.  They  can 
therefore  be  termed  random  gratings.  The  region  Xj  >  {(x, )  is  vacuum,  while  the  region 
*1  <  C(x,)  is  either  a  semi-infinite  metal  or  dielectric,  or  a  metal  or  dielectric  film  that 
is  either  free-standing  or  deposited  on  a  substrate.  The  surface  profile  function  C(x,)  is 
assumed  to  be  a  single-valued,  continuous,  differentiable  function  of  x , ,  and  to  constitute 
a  stationary,  Gaussian,  stochastic  process.  It  is  defined  by  the  properties  <C(x,)>  =  0 
and  <C(x,)((x',)>  =  i’W'dx,  -  x',|).  The  angle  brackets  here  denote  an  average  over 
the  ensemble  of  realizations  of  the  surface  profile.  The  form  of  the  surface  height 
correlation  function  ITflx,  -x", I)  will  be  specified  below  for  each  of  the  surfaces 
considered.  The  plane  of  incidence  (the  x,X3  plane)  is  perpendicular  to  the  generators 
of  these  gratings.  The  computational  method  used  in  the  calculation  of  the  angular 
distribution  of  the  intensity  of  the  incoherent  component  of  the  light  scattered  from,  or 
transmitted  through,  the  structures  studied  here  is  the  one  described  in  [4.7,12],  and 
extensions  thereof.  Thus,  the  random  grating  occupies  the  segment  (  -  f,/2,  L/2)  of  the 
X,  axis,  which  is  divided  into  N  -  L/Ax  equal  segments,  and  is  illuminated  from  the 
vacuum  side  by  a  Gaussian  beam  of  half-width  w.  The  half-width  ol  the  intercept  of 
this  beam  with  the  x,-axis  is  then  g  =  w/cos  60  where  60  is  the  angle  of  incidence.  The 
total  differential  reflection  or  transmission  coefficient  is  calculated  for  each  0!  different 
realizations  of  the  surface  profile  function,  and  the  results  are  averaged  over  this  number 
of  surfaces.  Finally,  the  coherent  contribution  is  subtracted  to  yield  the  incoherent 
contribution  to  the  mean  differential  reflection  or  transmission  coefficient. 


2.  Enhanced  backscattering  from  random  gratings  of  even  and  odd  symmetry 

If  enhanced  backscattering  from  moderately  rough,  reflecting  surfaces  is  due  primarily 
to  the  coherent  interference  of  each  doubly  reflected  optical  path  with  its  time-reversed 
partner,  we  should  expect  to  see  it  from  any  random  surface  that  can  multiply  scatter 
liglii.  The  earliest  computer  simulation  studies  of  enhanced  backscattering  [  4,7, 1 !  ]  were 
carried  out  for  reflecting  surfaces  with  large  rms  slopes  defined  by  a  surface  profile 
function  ((x,)  that  was  a  single-valued  function  of  x,  and  constituted  a  stationary. 
Gaussian  process,  characterized  by  a  Gaussian  surface  height  correlation  function 
•^(1*1  -  x',|)  =  <((x,){(x',)>/<{'(x,)>  =exp(-(x,  -x',)Va^),  where  a  is  the  trans¬ 
verse  correlation  length  of  the  roughness.  In  subsequent  work  it  was  shown  that  enhanced 
backscattering  is  obtained  from  reflecting  surfaces  with  large  rms  slopes  when  the  surface 
profile  function  ((x,)  is  no  longer  a  Gaussian-distributed  random  variable  [8],  It  has 
also  been  shown  that  enhanced  backscattering  occurs  for  forms  of  the  surface  height 
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correlation  functions  M'(|x,  —  x',|)  other  than  the  Gaussian  one  used  in  the  first  studies, 
when  C(X|)  is  a  Gaussian-distributed  random  variable  [14], 

Continuing  our  investigations  of  the  kinds  of  random  surfaces  that  give  rise  to 
enhanced  backscattering,  we  explore  here  the  consequences  of  relaxing  the  common 
assumption  that  the  surface  profile  function  is  a  stationary  stochastic  process.  We  do 
this  by  studying  the  scattering  of  p-  and  s-polarized  light  from  random  metallic  gratings 
with  large  rms  slopes  characterized  by  profile  functions  that  are  even  and  odd  functions 
of  X,.  The  study  of  the  enhanced  backscattering  of  light  from  random  surfaces 
characterized  by  even  profile  functions  also  allows  us  to  investigate  the  enhanced 
scattering  in  the  specular  direction  observed  in  the  angular  distribution  of  the  intensity 
of  the  incoherent  component  of  the  scattered  light  in  earlier  work  on  perfectly  conducting 
surfaces  by  Nieto-Vesperinas  and  Soto-Crespo  [9]. 

The  random  gratings  we  studied  were  constructed  in  the  following  way.  We  first 
constructed  a  random  surface  profile  function  C(x,).  obeying  Gaussian  statistics  defined 
by  the  properties  <((x,)>  =  0  and  <{(x,){(x',)>  =  5'exp(-(x,  -x',)*/a^),  by  the 
method  described  in  appendix  A  of  [12],  for  x,  in  the  interval  ( —  L/2,  L/2).  The  surface 
profile  functions  of  even  and  odd  symmetry  in  this  interval  were  then  constructed 
according  to  {,.„(x,)  =  [C(x,)  +  C(— x,)]/2,  respectively.  We  see  immediately  that  the 
surface  profiles  {,,„(x,)  defined  in  this  way  are  no  longer  stationary  random  processes, 
because  the  point  x,  =  0  is  a  distinguished  point. 

We  have  calculated  the  contribution  to  the  mean  differential  reflection  coefficient 
(drc)  from  the  incoherent  component  of  the  scattered  light,  as  a  function  of  the  scattering 
angle  for  light  of  both  p-  and  s-polarization  incident  on  a  silver  surface  whose  surface 
is  a  random  grating  of  even  and  odd  symmetry.  In  figure  1  we  present  our  results  for 
scattering  from  gratings  of  even  symmetry.  For  each  polarization  a  well  defined  peak 
is  observed  in  in  the  retro-reflection  direction,  $,  =  -20°.  In  addition, 

a  second  well  defined  peak  is  observed  in  the  specular  direction,  0,  =  20°.  This  is  the 
peak  we  call  the  specular  enhancement  peak.  We  emphasize  that  this  is  a  peak  in  the 
angular  dependence  of  the  intensity  of  the  incoherent  component  of  the  scattered 
light.  It  is  already  present  in  the  Kirchhoff  approximation,  i.e.  in  a  single-scattering 
approximation.  This  is  seen  in  the  results  presented  in  figure  2,  in  which  we  display  the 
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Figure  I.  The  incoherent  contribution  to  the  mean  one  for  the  scattering  of  a  beam  of  light 
frem  a  random  grating  of  even  symmetry  on  a  silver  surface.  Oq  -  20";  X  »  0.6127 ;<m; 
<((0)  *  -  17.2  +  i0.498i  6  «  t-4I42fim;  o  —  2  ftm;  L  «  25  6/jm;  g  =  6.4  ^m;  S  300; 
■  1000.  (fl)  p-polarization;  (6)  s-polarization. 
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Figure  Z  The  incoherent  contribution  to  the  mean  Dec  for  the  scattering  of  a  p-polarized 
beam  oflight  from  a  random  grating  of  even  symmetry  on  a  perfect  conductor.  The  parameters 
used  in  obtaining  these  results  are  those  usnl  in  obtaining  figure  I.  (o)  the  total  incoherent 
contribution  to  the  oact  (fi)  the  contribution  from  the  single-scallering  processes:  (c)  the 
contribution  from  the  pure  double-scattering  processes:  (d)  the  contribution  from  the 
single-  and  double-scattering  processes,  including  the  interference  terms. 


contribution  to  the  mean  diflerential  reflection  coefficient  of  p-polarized  light  incident 
on  the  same  surfaces  used  in  obtaining  figure  I  but  ruled  on  a  perfect  conductor,  for 
which  such  calculations  are  much  simpler  than  for  a  metal.  No  evidence  of  enhanced 
backscattering  is  seen  in  the  single  scattering  contribution  plotted  in  figure  2(h),  while 
a  well  defined  specular  enhancement  peak  is  present.  The  pure  double-scattering 
contribution  plotted  in  figure  2(c)  shows  an  enhanced  backscattering  peak,  as  well  as 
a  weak  specular  enhancement  peak.  These  results  are  consistent  with  the  picture 
of  enhanced  backscattering  as  a  multiple-scattering  phenomenon,  and  of  specular 
enhancement  as  a  predominantly  single-scattering  eflect.  In  figure  3  we  present  our 
results  for  the  scattering  of  light  of  p-  and  s-polarization  from  random  gratings  of  odd 
symmetry  on  a  silver  surface.  The  experimental  conditions  and  roughness  parameters 
are  the  same  in  this  case  as  were  assumed  in  obtaining  figures  I  and  2.  For  each 
polarization  a  well  defined  peak  is  observed  in  <dRp.,/50,>i„„t  in  the  retro-reflection 
direction,  6,  =  —20°,  but  there  is  no  peak  in  the  specular  direction,  6,  =  20°. 

These  results  can  be  understood  qualitatively  with  the  following  arguments.  If  we 
consider  only  single-scattering  contributions  in  the  scattering  of  light  from  an  arbitrary 
random  grating,  the  amplitude  components  interfering  in  the  specular  direction  arise 
from  the  so-called  specular  points.  These  are  points  on  the  surface  at  which  its  slope  is 
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Rsure  X  The  same  as  figure  I  but  for  a  random  grating  of  odd  symmetry. 


zero.  For  optically  rough  surfaces  the  random  heights  of  the  specular  points  have  large 
fluctuations  (compared  with  the  wavelength  X  of  the  incident  light),  and  the  relative 
phase  of  the  far-fleld  contributions  is  completely  random,  i.e.  uniformly  distributed  in 
the  interval  ( —  re,  n).  This  destroys  the  coherent  (or  specular)  component  of  the  scattered 
light.  The  situation  is  the  same  for  random  gratings  with  an  odd  profile.  However,  for 
random  gratings  with  an  even  profile,  there  are  pairs  of  contributions,  arising  from 


Figure  4.  An  experimenul  result  for  the  incoherent  contribution  to  the  mean  Die  for  the 
scattering  of  a  p<poIanzed  beam  of  light  from  a  random  grating  of  even  symmetry  on  a 
gold<oated  {d^,  »  500  A)  photoresist  film  of  mean  thickness  «  10  pm  deposited  on  a 
glass  substrate  of  thickness  *  5mm.  $9  «■  10*;  i  >  0.6328pm;  ••  1.64;  -  I.Sl; 

^  a  2.3  pm;  a  -  9.5  pro. 
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symmetric  specular  points  at  the  same  height,  at  x,  and  -x,,  which  interfere 
constructively  with  a  fixed,  non-random,  phase  difference.  This  increases  the  intensity 
of  the  scattered  light  in  the  specular  direction  by  a  factor  of  two  over  that  of  the 
background. 

Specular  enhancement  has  been  observed  experimentally.  In  figure  4  we  present  the 
mean  drc  for  the  scattering  of  p-polarized  light  from  a  random  grating  on  a  gold-coated 
photoresist  film  deposited  on  a  glass  substrate.  The  surface  is  sufficiently  rough  that  the 
coherent  component  of  the  scattered  light  is  negligible.  The  surface  consists  of  150 
different,  contiguous  segments,  each  of  which  is  200  pm  long,  and  each  of  which  is  an 
even  function  of  x,,  measured  from  its  midpoint.  The  width  of  the  incident  beam  is 
3  cm,  so  that  the  drc  plotted  is  equivalent  to  the  result  of  averaging  over  150  different 
realizations  of  a  random  surface.  A  well  defined  peak  is  observed  in  the  specular  direction, 
8,  =  10°.  However,  no  enhanced  backscattering  peak  is  observed  in  this  figure,  because 
the  roughness  parameters  are  such  that  the  overwhelming  contribution  to  the  drc  is 
due  to  single-scattering  processes,  i.e.  the  Kirchhoff  approximation  is  valid 


3.  Enhanced  backscattering  from  a  dielectric  film 

In  the  earliest  numerical  simulation  studies  of  the  scattering  of  p-polarized  light  from 
a  random  grating  on  a  semi-infinite  dielectric  medium  (BaSO^)  no  enhanced  back- 
scattering  was  observed  [II].  In  the  same  work  it  was  also  shown  that  if  the  index  of 
refraction  of  the  dielectric  medium  was  artificially  doubled  then  enhanced  backscattering 
was  observed  in  this  polarization  [II].  Subsequent  calculations  of  the  scattering  of 
s-polarized  light  from  the  same  surfaces  predicted  enhanced  backscattering  [  12].  It  was 
later  found  that  a  more  modest  increase  in  the  refractive  index,  namely  by  a  factor  of 
about  1.4- 1.5,  was  sufficient  to  induce  enhanced  backscattering  in  p-polarization  [12]. 

It  is  interesting,  then,  to  ask  whether  there  are  other  ways  to  induce  enhanced 
backscattering  in  the  scattering  of  p-polarized  light  from  a  random  grating  on  a  dielectric 
medium.  In  this  section  we  describe  a  way  of  doing  so,  which  relies  on  the  dielectric 
medium  being  in  the  form  of  a  thin  film  rather  than  semi-infinite. 

We  consider  the  case  in  which  a  film  of  a  nearly  transparent  dielectric  of  average 
thickness  d  is  deposited  on  a  perfectly  conducting  substrate.  The  interface  Xj  =  -d 
between  the  film  and  the  interface  is  assumed  to  be  planar.  The  interface  between  the 
film  and  the  vacuum  above  it  is  a  random  grating  defined  by  a  surface  profile  function 
{(x, ),  that  is  a  Gaussianly  distributed  random  variable  with  the  properties  <C(Jti  )>  =  0, 
<f(Xi)f(-Xi)>  =  i^exp(  —  (x,  -  x'l  )’/a').  This  structure  was  suggested  by  recent  work 
of  Jakeman  and  his  colleagues  who  showed  in  a  series  of  papers  that  light  scattered 
from  a  deep  random  phase  screen  placed  in  front  of  a  mirror  displays  a  strong  increase 
in  its  intensity  in  the  backscattering  direction  [21-23].  The  random  surface  in  our  case 
plays  the  role  of  the  deep  random  phase  screen,  while  the  perfect  conductor  is  the  mirror. 

We  have  carried  out  numerical  simulations  of  the  scattering  of  light  of  p-polarization 
from  this  structure  [13].  In  figure  5(a)  we  present  our  results  for  the  angular  distribution 
of  the  intensity  of  p-polarized  light  scattered  incoherently  from  a  dielectric  film  of  mean 
thickness  d  =  4.8  pm.  In  these  calculations  a  value  of  the  mean  film  thickness  d  was 
used  that  was  four  times  the  value  of  the  rhs  height  of  the  surface  roughness  6.  This 
was  because  for  such  large  values  of  the  film  thickness  the  probability  of  obtaining 
values  of  IC(x,)|  larger  than  d,  which  would  produce  'holes’  in  the  film,  is  reduced  to 
a  negligible  level.  A  sharp  peak  is  see.,  in  <dR,/d5,>i„„b  m  the  retro-reflection  direction 
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Figure  5.  iu)  The  incoherent  contribution  to  the  mean  drc  for  the  scattering  of  a  p-polanzed 
beam  of  iight  from  a  random  grating  on  the  surface  of  a  film  of  Ba50«  deposited 
on  a  planar,  perfectly  conducting  substrate.  Bq  5°;  ^  =  0.6328 pm.  =  1.628  -+  10.0003, 
d  =  4.8 pm:  S  *  t.2pm;  a  =  2pm;  L  =  25.6 pm;  g  “  6.4pm;  \  =  300;  *  1000.  (6)  The 

same  as  in  (ul  except  that  the  random  grating  is  on  a  semi-infinite  BaSO«  substrate. 


ff,  =  -5°.  For  comparison,  we  present  in  figure  5(b)  the  result  for 

obtained  for  exactly  the  same  experimental  conditions  and  roughness  parameters  except 

that  the  dielectric  is  now  semi-infinite.  No  enhanced  backscattering  is  seen  in  this  case. 

We  believe  that  the  enhanced  backscattering  we  observe  in  figure  5(a)  can  be 
explained  by  a  modification  of  an  argument  given  by  Jakeman  [21]  to  explain  the 
enhanced  backscattering  of  light  from  a  deep  random  phase  screen  placed  in  front  of 
a  plane  mirror.  When  the  perfect  conductor  is  placed  well  beyond  the  focusing  plane 
of  the  random  surface,  as  is  the  case  for  the  results  presented  here,  the  coherent  addition 
of  the  contribution  from  a  given  light  path  that  interacts  with  the  random  surface  at 
two  different  points  due  to  its  reflection  from  the  perfect  conductor,  and  its  time-reversed 
partner  leads  to  an  enhancement  of  up  to  a  factor  of  two  in  the  intensity  of  light  scattered 
into  the  backward  direction.  The  angular  width  of  this  peak  is  determined  by  a  transverse 
scattering  length,  which  is  the  characteristic  distance  between  the  points  on  the  random 
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Figarc  6.  An  experimental  result  for  the  incoherent  contribution  to  the  mean  drc  for  the 
scattering  of  a  p^polarized  beam  ^  U^t  from  a  random  grating  on  tbe  surface  of  a  dielectric 
film  deposited  on  a  planar  gold  substrate.  Oq  >  «  0.6328 ;im;  *  I.4I;  d  «  i.Sfon: 

S  !.08  /im;  a  ^  3.06 /tm. 

surface  through  which  the  light  paths  from  the  source  pass  and  repass  in  being  scattered 
into  the  retro-reflection  direction. 

The  effect  described  in  this  section  has  been  observed  experimentally.  In  figure  6  we 
present  experimental  results  for  the  mean  differential  reflection  coefI...ient  for  the 
scattering  of  p-polarized  light  from  a  dielectric  film  on  a  gold  substrate.  Although  gold 
is  not  a  perfect  conductor  under  the  conditions  of  the  experiment,  it  is  highly  reflecting, 
and  the  consequence  of  this  is  that  a  well  defined  peak  in  the  retro-reflection  direction, 
0,=  —  S°,  is  seen  in  the  angular  dependence  of  the  intensity  of  the  incoherent  component 
of  the  scattered  light.  Despite  the  differences  between  the  parameters  defining  the 
theoretical  and  experimental  results,  a  comparison  of  the  differential  reflection  coefficients 
presented  in  figures  5(a)  and  6  shows  them  to  be  qualitatively  very  similar,  and  even 
quantitatively  close. 


4.  Enhanced  transmission  through  thin  films 

Not  all  of  the  interesting  manifestations  of  weak  localization  in  the  interaction  of  light 
with  random  gratings  are  observed  in  reflection.  It  was  shown  recently  [24]  that  the 
angular  distribution  of  the  intensity  of  the  incoherent  component  of  p-polarized  light 
transmitted  through  a  thin  metal  film  surrounded  by  vacuum,  whose  illuminated  surface 
is  randomly  rough  while  the  back  surface  is  planar,  displays  a  well  defined  peak  in  the 
direction  of  transmission  that  is  directly  opposite  to  the  direction  of  specular  reflection 
of  the  incident  light  (the  antispecular  direction).  It  is  believed  that  'he  physical  origin 
of  this  effect  is  the  scattering,  by  the  surface  roughness,  of  the  surface  polaritons  in  the 


Enhanced  backacattering  and  tranamission  <4  light 


SI37 


film  excited,  through  the  roughness,  by  the  incident  light.  The  coherent  interference  of 
a  doubly  scattered  light/surface  polariton  path  with  its  time-reversed  panner  gives  the 
dominant  contribution  to  this  enhanced  transmission. 

The  original,  perturbation-theoretic  calculation  of  this  effect  was  extended  in  two 
directions  in  subsequent  numerical  simulation  calculations  [  14],  Rougher  surfaces  were 
studied,  and  calculations  were  carried  out  for  metal  films  on  semi-infinite  dielectric 
substrates.  In  addition,  the  phenomenon  of  enhanced  transmission  was  observed 
experimentally  [14], 

If  we  relax  the  assumption  that  only  the  illuminated  surface  of  the  film  is  a  random 
grating,  and  allow  the  back  surface  to  be  a  random  grating  as  well,  the  enhanced 
transmission  of  light  through  the  resulting  structure  can  be  amplified  with  respect  to 
what  it  is  when  only  the  illuminated  surface  is  rough.  In  figure  7(a)  we  present  results 
for  the  incoherent  contribution  to  the  mean  differential  transmission  coefficient  (dtc) 
for  the  transmission  of  p-polarized  light  through  a  thin  silver  film,  whose  illuminated 
surface  is  defined  by  xj  =  while  its  back  surface  is  given  by  Xy=—d  +  C(X|). 
Thus,  the  thickness  of  the  film  is  the  same  for  each  value  of  x,.  The  surface 
profile  function  ((x,)  is  a  stationary,  Gaussian  process  defined  by  <C(tC|)>  =  0  and 
<C(*i){(x',)>  =  i^exp(-(x,  —  x',)^/a^).  A  narrow,  high  peak  in  the  antispecular 
direction,  d,  =  -20°,  is  observed  in  the  angular  distribution  of  the  intensity  of  the 
incoherent  component  of  the  light  transmitted  through  the  film.  This  result  should  be 
compared  with  the  one  displayed  in  figure  7(h),  which  presents  the  corresponding  result 
for  a  film  whose  back  surface,  at  xj  =  —  d,  is  planar.  Although  the  overall  transmitted 
intensity  is  higher  in  this  case  than  in  the  case  of  a  film  with  both  surfaces  rough,  the 
enhanced  transmission  peak  is  broader  and  lower  compared  with  the  background  at 
its  position  that  it  is  in  the  latter  case. 

In  section  2  we  have  discussed  the  scattering  at  light  from  random  gratings  on  metal 
surfaces  when  the  surface  profile  function  of  the  random  grating  is  an  even  or 
odd  function  of  x,.  It  was  found  that  the  angular  dependence  of  the  intensity  of  the 


Rg««  7.  The  incoherent  contribution  to  the  mean  dtc  for  the  irifumiuion  of  a  p-polarized 
beam  of  light  through  a  thin  silver  him  with  rough  surfaces.  9^  «  20*.  i  »  0.4S79/im; 
((»)-  -'7,S°»'i0.24;  d-iSOOA,  ^-  282  8 A;  a  -  ISOOA;  L-  13.7^;  g  .  3.425/110; 
\  .  300;  1000.  (o)  The  illuminated  and  back  surfaces  are  at  Jts*C(Xi)  and 

X,  .  -d  >  ((x, ),  respectively;  (b)  the  illuminated  and  back  surfaces  are  at  Xg  .  C(X|)  and 
X  .  -d,  respectively. 
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incoherent  component  of  the  scattered  light  displayed  not  only  a  peak  in  the  retro- 
rehection  direction  but  also  a  peak  in  the  specular  direction  when  the  scattering  occurred 
from  random  gratings  characterized  by  even  surface  profiles.  Similar  effects  occur  in 
the  transmission  of  p-polarized  light  through  a  thin  metal  film  whose  illuminated  face 
is  a  random  grating  of  even  symmetry,  and  whose  back  face  is  planar.  In  figure  8  we 
plot  the  contnbution  to  the  mean  dtc  from  the  incoherent  component  of  the  transmitted 
light  when  light  of  p-polarization  is  incident  on  a  silver  film  whose  mean  thickness  is 
700  A,  and  whose  illuminated  surface  is  a  random  grating  of  even  symmetry.  From  figure  8 
we  see  that  the  transmission  of  p-polarized  light  thrpugh  a  metal  him  whose  illuminated 
surface  is  a  random  grating  of  even  symmetry  displays  a  peak  -i,  the  antispecular  direction 
6,  =  —do  (enhanced  transmission)  and  a  second  peak  in  the  forward  direction  8,  =  6o. 
which  we  describe  as  being  due  to  enhanced  refraction. 

Thus,  the  imposition  of  the  constraint  of  evenness  of  an  otherwise  random  surface, 
and  the  consequent  removal  of  its  stationarity,  does  not  eliminate  the  enhanced 
transmission  that  occurs  in  the  transmission  of  p-polarized  light  through  a  metal  him 
with  a  random  surface  that  has  no  such  symmetry  property,  and  gives  rise  to  a  new 
effect— enhanced  refraction.  The  phenomenon  of  enhanced  refraction  is  believed  to  be 
due  to  the  coherent  interference  of  light  scattered  from  correlated  pairs  of  points  on  the 
random  surface,  which  already  occurs  in  a  single-scattering  approximation. 

Enhanced  refraction  has  now  been  observed  experimentally.  In  hgure  9  we  present 
the  angular  distribution  of  the  intensity  of  the  incoherent  component  of  p-polarized 
light  transmitted  through  a  gold-coated  photoresist  plate  of  mean  thickness  d  =  10  pm 
deposited  on  a  glass  him  of  5  mm  thickness.  It  was  fabricated  in  the  manner  described 


Figure  a.  The  incoherent  coninbulion  to  the  mean  ore  tor  the  transmiuion  of  a  p-polarized 
beam  of  tight  through  a  silver  film  whose  illuminated  face  is  a  random  grating  of  even 
symmetrv  while  its  back  face  ts  planar  Fo“5'.  /.»  0.4579 ^m;  elcu)  -  — 7.5  +  i0.24; 
4  •  700  A.  li-  282  84  A.  u  ^  ISOOA;  I  -  I3  7um;p  •  3.425 pm;  N  -XO.S,-  1500 
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Figure  9.  An  enpenmental  result  for  the  incoherent  contribution  to  the  mean  dtc  for  the 
transmission  of  a  p-polarized  beam  of  light  through  a  random  grating  of  even  symmetry 
on  a  gold-coated  (d^,  •  SOOAl  photoresist  Aim  of  mean  thickness  d^  »  lOitm  deposited 
on  a  glass  substrate  of  thickness  tf^  =  5  mm.  Aq  *  20\/  *  0.6328  pm;  =  1.64.  n,,  =  1.51. 
S  ■  1,4  pm.  a  *  8.0  pm 


in  section  2.  A  prominent  peak  is  seen  in  the  angular  dependence  of  the  transmitted 
intensity  in  the  forward  direction,  i.e.  for  a  scattering  angle  S,  =  10°.  It  is  due  to  enhanced 
refraction.  The  absence  of  enhanced  transmission  in  this  case  is  consistent  with  results 
of  numerical  simulations  which  show  that  enhanced  transmission  is  suppressed  as  d/a 
is  increased,  presumably  because  the  surface  polaritons  responsible  for  the  effect  [24] 
are  strongly  damped  by  the  roughness.  The  subsidiary  maxima  in  the  mean  dtc  are 
believed  to  be  due  to  multiple  reflections  from  the  two  interfaces  in  the  structure  studied. 


S.  Conclusions 

The  results  presented  in  this  paper,  taken  together  with  those  of  earlier  investigations 
[8, 1 4],  show  that  the  enhanced  backscattering  of  light  from  moderately  rough,  reflecting, 
random  surfaces  that  scatter  light  multiply  is  a  remarkably  robust  phenomenon.  Provided 
that  these  surfaces  scatter  light  multiply,  it  is  independent  of  the  statistics  of  the  surface 
profile  function;  it  is  independent  of  the  form  of  the  surface  height  correlation  function; 
and  it  occurs  for  surfaces  for  which  the  surface  profile  function  is  not  a  stationary 
stochastic  process.  These  results  strongly  suggest  that  the  dominant  property  of  a 
moderately  rough,  reflecting  surface  that  is  responsible  for  enhanced  backscattering  is 
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its  ability  to  scatter  the  incident  light  multiply:  other  properties  of  such  surfaces  may 
affect  the  form  the  differential  reflection  coefficient  possesses,  but  not  the  existence  of 
enhanced  backscattering  itself. 

In  addition,  when  the  random  surface  profile  is  constrained  to  be  an  even  function 
of  X,,  a  new  scattering  effect  is  observed,  specular  enhancement,  which  is  absent  in 
general  in  the  scattering  from  an  unconstrained  random  surface  (see,  however,  [8]). 
This  effect  however,  is  already  present  in  a  single-scattering  approximation,  and  thus  is 
not  primarily  a  multiple-scattering  effect. 

The  results  presented  in  section  3  show  that  a  dielectric  medium  bounded 
by  a  moderately  rough  random  grating  surface  which  does  not  display  enhanced 
backscattering  of  p-polarized  light  when  it  is  semi-infinite  can  be  made  to  do  so  by 
depositing  the  dielectric  medium  in  the  form  of  a  thin  film,  whose  interface  with  vacuum 
is  a  random  grating  of  large  rms  slope,  on  the  planar  surface  of  a  perfect  conductor. 
The  observed  enhanced  backscattering  from  this  structure  is  believed  to  be  due  to  the 
coherent  interference  of  a  light  path  that  passes  through  the  random  surface  twice,  due 
to  its  reflection  from  the  perfect  conductor,  and  its  time-reversed  partner.  In  effect,  the 
presence  of  the  perfect  conductor  ensures  the  douole-scattering  of  the  incident  light 
from  the  random  surface  that  is  primarily  responsible  for  enhanced  backscattering. 

Finally,  the  enhanced  transmission  of  p-polariaed  light  through  a  thin  metal  film, 
whose  back  surface,  as  well  as  its  illuminated  surface,  is  a  random  grating  is  noticeably 
stronger  than  it  is  when  the  back  face  is  planar.  This  conclusion  may  make  this  effect 
easier  to  observe  experimentally.  When  the  back  face  of  such  a  film  is  planar,  and  the 
illuminated  face  is  a  random  grating  characterized  by  a  surface  profile  function  that  is 
an  even  function  of  x,,  the  incoherent  component  of  the  transmitted  light  displays 
enhanced  refraction  in  addition  to  enhanced  transmission.  The  former  effect,  however, 
IS  pnmarily  a  single-scattering  effect,  in  contrast  with  the  multiple-scattering  nature  of 
the  latter. 
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Abstract  Measurements  are  presented  of  the  angular  distribution  of  four  wavelengths  of 
light  scattered  by  a  one>diinensional  random  rough  surface,  whose  probability  density 
function  is  Gaussian  with  a  standard  deviation  a  -  \22  ±  0.02  pm  and  whose  lateral 
correlation  function  is  also  Gaussian  with  1  je  width  t  «  3  1 7  +  0  07  pm  The  wavelengths 
used  are  063.  115.  3.39  and  t06pm.  The  surface  is  used  in  two  forms  coated  with  gold 
and  as  an  almost  lossless  dielectric.  The  results  are  compared  to  those  predicted  by  a 
double  scattering  form  of  the  KirchholT  formulation.  Agreement  is  good  at  small  angles  of 
incidence  but  less  good  at  larger  angles  of  incidence. 


1.  Introduction 

The  experimental  observation  of  enhanced  backscattering  from  random  rough  surfaces 
of  large  root-mean-square  slope,  first  reported  by  Mendez  and  O'Donnell  [1,2],  has 
stimulated  a  re-examination  of  the  problem  of  light  scattering  in  the  past  few  years.  The 
main  progress  to  date  has  been  the  development  of  ‘exact'  numerical  codes  for  the 
solution  of  Maxwell’s  equations  from  a  one-dimensional  surface  illuminated  with  either  s 
(i.e.  TE)  polarization  or  p  (i.e.  tm)  polarization  [3-7].  With  s  polarization,  the  electric 
vector  is  parallel  to  the  grooves,  whereas  with  p  polarization  it  is  perpendicular  to  the 
grooves,  as  in  figure  I  (this  figure  also  shows  the  sign  convention  used  for  the  incident 
and  scattering  angles).  An  important  feature  of  the  work  of  Mendez  and  O'Donnell 
was  that  the  surfaces  were  relatively  well  characterized,  with  Gaussian  statistics  for  the 
surface  height  and  a  single-scale  Gaussian  correlation  function.  Since  the  statistics  of 
the  surface  were  known,  a  critical  comparison  between  experiment  and  theory  could  be 
made  with  confidence. 

The  shape  of  the  scattering  cross  section  curves  with  angle  of  observation  for 
high-sloped  surfaces  is  quite  different  from  that  for  simple  low-sloped  ones  and  early 
numerical  results  [3]  were  encouraging  since  they  were  in  fairly  good  agreement  with 
the  experimental  ones  particularly  at  near-normal  incidence.  In  order  to  carry  out  a 
more  critical  comparison  between  real  experiments  and  numerical  ones,  it  is  important 
that  the  surface  is  well  characterized  and  also  helpful  if  a  range  of  wavelengths  are  used. 
The  results  presented  here  are  intended  to  supplement  those  already  reported  [8-10] 
with  the  aim  of  providing  a  reliable  body  of  experimental  data  for  comparison  with 
numerical  work  and  analytical  theory.  The  surface  used  is  one-dimensional,  for  two 
reasons:  firstly,  it  can  be  characterized  much  more  accurately  than  a  two-dimensional 
one,  since  a  sharp  chisel-shaped  stylus  can  be  used  in  a  surface  profilometer;  secondly. 
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Figure  I.  Polarization  and  angle  notation  Tor  in-piane  scattering  from  a  one-dimcnsional 
rough  surface. 


‘exact’  numerical  calculations  of  the  light  scattering  are  only  feasible  at  the  present  time 
for  the  one-dimensional  case. 

When  comparing  experimental  measurements  of  light  scattering  with  numerical 
computations,  it  is  helpful  if  the  numerical  results  provide  some  physical  insight  to  the 
scattering  process.  For  example,  it  is  believed  from  the  experiments  that  the  mechanism 
giving  rise  to  the  enhanced  backscatter  peak  is  multiple  scattering;  however,  numerical 
calculations  based  on  exact  theory  do  not  separate  the  single  and  multiple  scatter  terms, 
and  therefore  do  not  provide  the  insight  that  is  desirable  (however,  iterative  solutions 
do  separate  the  single  and  multiple  scatter  terms).  For  this  reason,  we  have  written 
numerical  code  based  on  a  multiple  (double  and  triple)  scattering  extension  of  the 
Kirchhoff  boundary  condition,  including  the  effects  of  shadowing  (see  [  1 1  ]  for  details 
and  further  references).  In  section  3  of  this  paper  we  compare  the  results  of  this  code 
with  the  experimental  results  and  ‘exact’  numerical  code. 


2.  Experimental  results 

Master  surfaces  are  produced  by  exposing  a  thick  layer  of  photoresist  (%  12 /jm  of 
Shipley  S1400-37)to  several  statistically  independent  laser  speckle  patterns.  Two  versions 
of  the  surface  were  prepared  using  a  replication  technique  [8],  one  being  coated  with 
%  1000  A  of  gold  and  the  other  being  an  almost  lossless  dielectric  of  refractive  index 
n  =  1.41  (at  X  =  0.63  pm).  Figure  2  shows  the  probability  histogram  of  surface  height 
and  surface  autocorrelation  function,  as  measured  by  a  Talystep  prohlometer  whose 
stylus  is  a  pyramid  of  70°  apex  angle  truncated  by  a  flat  region  of  %0.S  pm.  Both  are 
good  fits  to  Gaussian  functions,  with  the  root-mean-square  height  tr  =  1.22  ±  0.02  pm 
and  1  /e  correlation  length  r  =  3.17  +  0.07  pm.  The  angular  distribution  of  the  scattered 
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Figure  1  Probability  histogram  (upper)  and  autocorrelation  function  I  lower  i  of  the  surface 
height  fluctuation  as  calculated  from  Talysiep  measurements  for  the  surface  i  #  46 1  used  in 
this  paper. 


light  was  measured  at  four  wavelengths  (0.63 grm,  1.15  frm,  3.39  fim  and  106  pm)  using 
the  equipment  described  in  [  10].  For  each  angle  of  incidence,  measurements  are  made 
with  p-polarization  incident  and  p-polarization  collected  Cp-p’  scattering)  and  s- 
polarization  incident  and  s-polarization  collected  (‘s-s’  scattering);  no  crossed  polarized 
components  were  detectable.  For  a  perfect  conductor,  these  measurements  give  a 
complete  description  of  the  scattering  characteristics  of  the  surface,  but  in  general  four 
scattering  coefficients  are  required  for  materials  of  finite  conductivity;  these  can  be  found 
by  measuring  the  polarization  of  the  scattered  light  for  various  input  polarizations.  Also, 
the  measurements  reported  here  yield  the  relative  scattering  cross  section,  as  no  absolute 
calibration  is  performed. 

The  relative  scattering  cross  sections  for  angles  of  incidence  of  0°,  -30°  and  -60° 
and  the  four  wavelengths  are  shown  in  figures  3  and  4  for  the  gold-coated  surface  and 
figures  5  and  6  for  the  dielectric  surface.  The  enhanced  backscatter  peak,  where  present, 
occurs  on  the  right-hand  side  of  the  graphs  (i.e.  at  positive  angles,  see  figure  1  for  the 
sign  convention  for  the  angles)  and  any  specular  component  is  on  the  left-hand  side 
(i.e.  negative  angles);  for  the  10.6pm  measurements,  the  specular  peak  was  very  much 
greater  than  the  diffuse  component  and  is  not  shown.  A  few  features  are  of  particular  note. 
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Figure  3.  Relative  scattering  cross  section  as  a  function  of  scattering  angle  for  the  gold-coated 
surface,  for  angles  of  incidence  of  O',  -30'  and  -60®.  for  p-p  scatienng  (open  circles)  and 
5-s  scattering  (crosses).  The  iefl-band  column  is  for  a  wavelength  of  0.63 for  which 
<j/x  as  1.93  and  t/A  s  5.02  and  (he  right-hand  column  for  A  =  1.15/im,  for  which  <r/A  =  J.07 
and  t/A  =  2.76.  The  enhanced  backscatter  peak,  where  present,  occurs  at  positive  angles 
(right-hand  side  of  each  graph). 


( 1 )  The  enhanced  backscatter  peak  and  sidelobe  structure  are  clearly  visible  for  the 
shorter  wavelengths  at  an  angle  of  incidence  less  than  approximately  —30°  for  the 
gold-coated  surface;  the  width  of  the  peak  is  proportional  to  the  wavelength.  The  peak 
is  not  observed  for  the  p-p  scattering  at  10.6 pm  for  the  gold  surface  or  for  scattering 
from  the  dielectric. 

(2)  With  the  exception  of  the  p-p  case  at  10.6  pm,  the  p-p  and  s-s  scattering  by  the 
gold  surface  are  very  similar;  for  the  dielectric  surface,  however,  the  p-p  and  s-s  scattering 
cross  sections  are  quite  different,  as  one  might  e-pect  by  analogy  with  reflection  from 
a  planar  surface.  Using  a  value  of  n  =  1.41  for  the  refractive  index  of  the  (almost  lossless) 
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Figure  4.  As  for  figure  3  but  wavelengths  of  3.39 ;im  (o,/  *  0.36,  ilk  -  094)  and  l0.6/im 
(cr/x  »  0,12,  t/A  »  0.30).  For  the  10.6 /jm  curves,  the  angles  of  incidence  were  0  .  -  20  and 
-  40®.  The  ( strong)  specular  component  in  the  106  curves  is  not  shown 


dielectric  gives  a  Brewster  angle  of  5;  55°.  Considering  single  scattering  to  be  the  dominant 
mechanism  and  treating  this  as  a  reflection  from  a  locally  plane  surface  gives  an 
expected  minimum  of  the  p-p  scattered  intensity  at  an  angle  equal  to  approximately 
(-110°-  incident  angle):  the  angles  are  roughly  in  accordance  with  this  simple  picture. 
The  s-s  and  p-p  scattered  intensities  in  the  backscatter  direction  appear  to  be  almost 
equal  to  each  other  for  all  angles  of  incidence  and  wavelengths,  for  the  dielectric. 

(3)  The  overall  shape  of  the  curves  is  dramatically  different  from  the  Gaussian-type 
shapes  (centred  on  the  specular  angle)  normally  encountered  in  scattering  from 
low-sloped  surfaces. 

The  principal  purpose  of  figures  3  to  6  is  to  provide  a  reliable  set  of  data  for 
comparison  with  numerical  calculations,  and  analytical  theories  should  any  become 
available. 
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Figure  5.  As  for  figure  3  but  for  the  dielectric  surface.  *  1.4 1,  n,  ,j  *  1  40  The  vertical 
scale  is  not  the  same  as  that  used  in  figure  3  (both  are  relatii'e  scattering  cross  sections). 


3.  KirchholT  multiple  scatter  approximation 

One  can  compare  the  above  experimental  results  to  those  of'exact’  numerical  calculations 
based  on  the  extinction  theorem  and  its  extensions  [3-7],  and  some  comparisons  of 
experiment  and  calculations  for  a  perfect  conductor  were  given  in  [10].  Although  such 
comparisons  are  valuable,  one  problem  with  the  ‘exact’  numerical  solution  is  that  it 
gives  little  physical  insight  into  the  problem.  Wc  have  therefore  attempted  to  extend 
the  Kirchhoff  approximation  (i.e.  tangent  plane  approximation  for  each  scattering  event) 
to  double  (and  multiple)  scattering  [11]. 

The  numerical  calculations  were  carried  out  using  the  method  descrit>ed  in  [  1 1  ]  for 
a  perfect  conductor;  typically  the  energy  conservation  (unitarity)  held  to  better  than 
3%  considering  just  the  single  and  double  scatter  terms  for  surface  #  46  (except  for  the 
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Figure  6.  As  for  figure  4  but  for  (he  dielectnc  surface,  3^  -  *^10  6  ^ 


—  60°  results  for  which  the  departure  from  unitarity  is  5: 6% )  and  the  results  are  averaged 
over  approximately  10^  realizations.  Figures  7  and  8  show  the  results  of  the  calculations 
for  /.  =  0.63  ;rm  and  1 . 1 5  /rm  respectively,  for  incident  angles  of  0°,  -  30°  and  -  60°  and 
s-s  and  p-p  scattering.  Each  graph  shows  the  single,  double  and  total  scattered  intensity. 
The  enhanced  backscatter  peak  occurs  only  in  the  double  scatter  component,  showing 
conclusively  that  the  enhancement  is  a  multiple  scattering  effect.  The  enhancement  is 
on  the  order  of  a  factor  of  two  in  the  double  scattered  component  for  all  angles  of 
incidence,  but  the  enhancement  in  the  total  intensity  is  much  less  than  two  and  decreases 
with  increasing  angle  of  incidence  due  to  the  fact  that  the  double  scattered  intensity 
also  decreases  with  incidence  angle. 

Figure  9  compares  the  total  scattered  intensity  for  s-s  scattering  from  figures  7  and 
8  with  the  results  of  'exact'  numerical  calculations  (based  on  the  extinction  theorem 
method  for  a  perfect  conductor  [3])  and  the  experimental  results  of  figure  3,  for 
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Figure  7.  Numerically  calculated  scattered  intensities  for  a  perfect  conductor,  using  the 
doubIC'Seattering  Kirchholf  approximation,  for  angles  of  incidence  of  O'.  —30'  and  —60 
and  S'S  scattenng  (leftj  and  p-p  scattering  (nght).  TTie  wavelength  is  0.63  ^m.  Each  graph 
shows  the  doubly  scattered  intensity  (lowest  curve),  single  scattered  (middle)  and  total 
intensity  (coherent  sum)  (top  curve). 


/  =  0.63  >im  and  1.15  at  three  angles  of  incidence.  The  two  numerical  calculations 
agree  well,  showing  that  the  KirchhofT  approximation  is  reasonable  for  these  surface 
parameters  (the  average  radius  of  curvature,  defined  as  the  inverse  of  the  standard 
deviation  of  the  surface  curvature  2^0 Ix^,  is  »  2.4  pm  for  surface  #  46)  and  both 
agree  well  with  the  experimental  measurements  for  zero  angle  of  incidence.  However, 
there  is  a  clear  discrepancy  between  experiment  and  numerical  calculation  for  the  -  30° 
and  -60°  angles  of  incidence.  (This  genera)  behaviour  is  also  shown  in  the  case  of  p-p 
scattering.) 
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Figure  8.  As  for  figure  7,  but  wavelength  of  1.15/im. 


One  possible  cause  of  the  discrepancy  could  be  that  the  calculations  are  for  a  perfect 
conductor,  whereas  the  experiments  are  for  a  real  metal  (gold).  However,  calc  ilations 
by  ourselves  and  others  [S]  show  that,  for  these  values  of  rms  surface  height  and 
correlation  length,  there  is  very  little  difference  between  the  results  for  gold  and  for  a 
perfect  conductor,  particularly  for  s-s  scattering.  One  problem  with  most  methods  of 
calculation,  including  that  used  here,  is  that  a  very  small  length  of  surface  is  considered, 
giving  rise  to  the  possibility  of  an  ‘end-effect’  error  (e.g.  due  to  long-range  surface 
plasmons);  however,  the  method  of  calculation  of  Saillard  and  Maystre  [7]  uses  an 
extremely  long  length  of  surface  with  good  agreement  with  the  other  calculations  and 
poor  agreement  with  the  measurements  at  larger  angles  of  incidence. 

It  seems,  therefore,  that  there  may  be  some  aspect  of  the  experiment  that  does  not 
correspond  to  the  calculations.  Previous  results  for  a  Lambertian  diffuser  have 
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Figure  9.  Comparison  of  Kirchhoff  calculation  (from  figures  7  and  8).  'exact'  numerical 
calculation  and  experiment  (from  figure  3).  for  5*5  scattering  at  a  «  and  l.lS/im, 

and  angles  of  incidence  equal  to  0\  -30'’  and  —60’.  The  solid  curves  arc  the  Kirchhoff 
calculation,  triangles  the  'exact'  calculation  and  crosses  are  the  expenmenial  results.  Note 
the  good  agreement  between  the  two  numerical  calculations  but  the  departure  of  the 
experimental  results  for  larger  angles  of  incidence 


demonstrated  that  the  scalterometer  measures  the  correct  quantity  [10].  The  measure- 
mem  of  surface  properties  might  be  in  error.  If  one  calculates  the  scattered  intensity 
for,  say,  -60°  angle  of  incidence  for  a  surface  that  has  an  rms  roughness  50%  larger 
than  the  measured  value,  then  reasonable  agreement  is  obtained  between  experiment 
and  numerical  calculation.  However,  (a)  it  is  extremely  unlikely  that  such  a  gross  error 
could  occur  (stylus  tips  effects  a.,  discussed  by  Church  [12])  and  (fc)  the  agreement 
for  0'  angle  of  incidence  is  then  very  poor  indeed,  particularly  as  regards  the  location 
of  the  minima  around  the  backscatter  peak.  Ishimaru  and  Chen  [13]  have  shown  that 
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a  departure  from  Gaussianit>  of  the  correlation  function  could  be  responsible  for  the 
discrepancy,  and  the  measured  correlation  does  show  a  small  departure  from  the 
Gaussian  shape.  However,  it  is  notoriously  difhcull  to  estimate  the  correlation  function 
of  stylus  traces  and  the  departure  shown  in  figure  2  is  characteristic  of  inadequate 
de-trending  of  the  mean;  the  method  of  manufacture  of  the  surfaces  strongly  encourages 
a  Gaussian  correlation  of  surface  height.  The  cause  of  this  discrepancy  for  larger  angles 
of  incidence  is  therefore  not  resolved  at  the  present  time. 


4.  Summary 

A  set  of  scattering  data  for  a  one-dimensional  surface  at  four  wavelengths,  three  angles 
of  incidence  and  two  materials  has  been  presented  for  critical  comparison  with  numerical 
calculations  and  theoretical  studies.  A  multiple  scatter  extension  of  the  KirchholT 
approximation  has  been  shown  to  provide  additional  physical  evidence  that  the 
predominant  cause  of  the  enhanced  backscatter  peak  is  due  to  multiple  scattering.  There 
remains  a  significant  disagreement  between  experiment  and  numerical  calculations  for 
large  angles  of  incidence  the  cause  of  which  is  still  unresolved 

The  data  presented  in  figures  3  to  6,  together  with  sample  Talystep  traces,  is  available 
on  a  PC-  or  Macintosh-compatible  diskette  on  application  to  the  first  author 
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